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FOREWORD 


The Quarterly of Applied Mathematics has been founded primarily to 
meet the needs of certain mathematicians and engineers whose interests ex- 
tend beyond the accepted boundaries of their respective groups. These mathe- 
maticians find their greatest interest in the application of mathematics to 
physical problems, and these engineers seek solutions of practical problems 
by advanced mathematical methods. Thus they meet on the common ground 
of applied mathematics with a stimulating variety of interest. 

It is not desirable to attempt too »recise a definition of the boundaries of 
the field to which the Quarterly wil’ devoted. The mathematical solution 
of one problem often throws light on another problem in an entirely different 
field; indeed, the peculiar strength of the mathematical method lies in its 
power to cut across those lines of demarcation which seem to divide science 
into separate compartments. 

Nevertheless, it is necessary to giv’ an outline of policy for, within fairly 
wide limits, the pages of the Quarterly .aould appeal to a common interest. It 
seems best to start with the common ground of mathematics and engineering as 
a nucleus, and to build around it a wider circle of interest, embracing mathe- 
matical theory related to engineering problems. Thus certain subjects—fluid 
mechanics, elasticity, plasticity, thermodynamics, and classical mechanics in 
its engineering applications—are to be regarded as lying within the scope of 
the Quarterly, and to these must be added electrical engineering, which has 
been one of the most fruitful fields of mathematical application. 

While it is not the purpose of the Quarterly to publish experimental re- 
sults, we shall welcome mathematical contributions which have an intimate 
connection with application in industry or practical science. Indeed, the ideal 
contribution to our pages would be one in which advanced and general mathe- 
matical methods lead speedily to results which are in close agreement with 
experiment, and which are of high importance, either in direct practical ap- 
plication or as an illumination of interesting phenomena hitherto unex- 
plained. 

THE EpITors. 








TOOLING UP MATHEMATICS FOR ENGINEERING* 


BY 
THEODORE von KARMAN 


California Institute of Technology 


It has often been said that one of the primary objectives of Mathematics 
is to furnish tools to physicists and engineers for solution of their problems. 
It is evident from the history of the mathematical sciences that many funda- 
mental mathematical discoveries have been initiated by the urge for under- 
standing nature’s laws and many mathematical methods have been invented 
by men primarily interested in practical applications. However, every true 
mathematician will feel that a restriction of mathematical research to prob- 
lems which have immediate applications would be unfair to the “Queen of 
Sciences.” As a matter of fact, the devoted “minnesingers” of the Queen have 
often revolted against degradation of their mistress to the position of a “hand- 
maiden” of her more practical minded and temporarily more prosperous 
sisters. 

It is not difficult to understand the reasons for the controversial view- 
points of mathematicians and engineers. They have been pointed out more 
than once, by representatives of both professions. 

The mathematician says to the engineer: | have built a building on a sound 
foundation: a system of theorems based on well defined postulates. I have 
delved into the analysis of the process of logical thinking to find out whether 
or not there are any statements which could be considered true or at least 
potentially true. I am interested in functional relations between entities 
which are well defined creations of my own mind and in methods which en- 
able me to explore various aspects of such functional relations. If you find 
any of the concepts, logical processes or methods which I have developed 
useful for your daily work, I am certainly glad. All my results are at your 
disposal, but let me pursue my own objectives in my own way. 

Says the engineer: Your great forbears, who were mathematicians long 
before you, talked a different language. Did not Leonhard Euler distribute 
his time between discoveries in pure mathematics and in the theory of engi- 
neering devices? The fundamentals of the theory of turbines, the theory of 
buckling of columns, the theory of driving piles into soil were contributions 
of Euler. The development of mathematical analysis cannot be separated 
from the development of physics and especially of mechanics. It is doubtful 
whether a human mind would ever have conceived the idea of differential 
equations without the urge to find a mathematical tool for the computation 
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of the path of moving bodies. If one assumes that the motion is determined by 
certain fundamental mechanical or geometrical relations, which are valid at 
every instant of the motion, one naturally is led to the idea of the differential 
equation. Also, the calculus of variations was invented mainly for solution 
of physical problems; some of which were of teleological, some of practical 
nature. The eighteenth century and the first decades of the nineteenth were 
perhaps the period of the most glorious progress in mathematical science; at 
that time, there was no distinction between pure and applied mathematicians. 
The abstract minded mathematicians stepped in after the big job was done; 
they endeavored to fill certain logical gaps, to systematize and codify the 
abundance of methods and theorems which the giants of the foregoing period 
created by a combination of logical thinking and creative intuition. 

The mathematician: It seems to me that you underestimate the impor- 
tance of what you call systematization and codification. Don’t you think that 
in order to assure the correct application of calculus and differential equa- 
tions, there was an absolute necessity to define exactly what we mean by a 
limiting process; or, was it not absolutely necessary to give a real sense to 
such terms as infinitely small and infinitely large? You may remember that 
Galileo—whom you hardly can call an abstract or pure mathematician— 
pointed out the contradictions which are unavoidable if you try to apply the 
notions of equality and inequality to infinite quantities. He noticed that you 
can say either that the number of the integers is larger than the number of 
the squares, since every square is an integer, but not every integer is a square; 
or you can say with the same justification that there are as many squares as 
integers, since every number has a square. The notions of commensurability, 
denumerability, the logical analysis of the continuum, the theory of sets, and 
in more recent times, topology, were fundamental steps in the development 
of the human mind. Many of these developments were conceived independ- 
ently of any conscious physical applications. But even for the sake of ap- 
plications, it was necessary to improve the foundations of our own house, 
that is to improve the logical structure of mathematics. Without exact analy- 
sis of the conditions for the convergence of series (the conditions which 
allow carrying out the processes of differentiation and integration), nobody 
could feel safe in handling series. It is not correct that the tendency to seek 
for a solid foundation of the new discoveries began after the men endowed 
with imagination and intuition did the big job. D’Alembert already de- 
manded that the calculus be founded on the methods of limits. Cauchy, 
Legendre and Gauss certainly were among the creative mathematical ge- 
niuses in your sense; they effectively contributed to the transition from intui- 
tion to rigor. In the second half of the nineteenth century this development 
continued toward the great goal that the mathematicians of that age—per- 
haps optimistically—considered as perfect logic or absolute rigor. However, 
in addition to the clarification of the fundamentals, that period also opened 
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new paths for applied mathematics. You mentioned, for example, differential 
equations. Don’t you believe that the theory of functions of complex vari- 
ables, the classification of differential equations according to their singu- 
larities, and the investigation of these singularities, all developed in the period 
that you call the period of codification, were most important steps in building 
up the very branch of mathematics from which you engineers derive so much 
benefit? These theories changed the primitive way of finding solutions of 
differential equations by trial into a systematic method of mastering the 
whole field. 

The engineer: I agree, especially with what you say about complex vari- 
ables. Indeed, the conformal transformation is one of the most powerful and 
most elegant methods for the solution of innumerable physical and engineer- 
ing problems. I also agree with you on the fundamental importance of the 
analysis of singularities. In fact, our graphical and numerical methods neces- 
sarily fail or become awkward near the irregular points and we have to take 
recourse to analytical methods. However, you mathematicians unfortunately 
are somewhat like a physician who is less interested in the laws of normal 
functioning ot the human body than in its diseases, or like the psychologist 
who instead of investigating the laws of normal mental! processes concen- 
trates his attention on the pathological aberrations of the human mind. We 
have to deal in most cases with “sound functions” and would like to have 
efficient methods to determine with fair accuracy their behavior in certain 
definite cases. 

Answers the mathematician: Can you not apply the general methods that 
we developed for the solution of differential and integral equations? If the 
solutions are given by “sound functions,” as you please to call them, I do not 
see any great difficulty nor do I see what more you expect us to do. 

The engineer: Your general theorems deal mostly with the existence of 
solutions and the convergence of your methods of solution. You may recall 
the wisecrack of Heaviside: “According to the mathematicians this series is 
divergent; therefore, we may be able to do something useful with it.” You 
people spend much time and much wit to show the existence of solutions 
whose existence often is evident to us for obvious physical reasons. You 
seldom take the pains to find and discuss the actual solutions. If you do so, 
then you restrict yourself mostly to simple cases, as for example, problems 
involving bodies of simple geometrical shapes. I refer to the so-called special 
functions. I concede that a great many such functions were investigated by 
mathematicians. Their values have been tabulated, their developments in 
series and their representations by definite integrals have been worked out 
in great detail. Unfortunately, such functions have only a restricted field of 
application in engineering. The physicist in his search for fundamental laws 
may choose specimens of simple geometrical shapes for his experimentation. 
The engineer has to deal directly with structures of complicated shapes; he 
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cannot give to a structure a simple geometrical form just because the stress 
distribution in such a structure can be calculated by special functions. Fur- 
thermore, most special functions are applicable only to linear problems. In 
the past, physicists and engineers often linearized their problems for simplic- 
ity’s sake. Mathematicians liked this simplification because it furnished a 
beautiful hunting ground for the application of elegant mathematical meth- 
ods. Unfortunately, as engineering science progressed, the need for more 
exact information and the necessity to get nearer and nearer to physical 
reality, forces us to grapple with many nonlinear problems. 

The mathematician: Well, many modern mathematicians are extremely 
interested in non-linear problems. It seems your primary need is the develop- 
ment of appropriate methods of approximation. However, you are not right 
in your criticism of our proofs of existence. Many proofs of existence in mod- 
ern mathematics go far beyond the limits of intuition. Then, too, I under- 
stand you engineers have good success with various iteration methods. Now, 
if we want to prove for example the existence of a solution of a boundary 
value problem, very often we use the iteration method. In other words, we 
really construct a sequence of approximate solutions exactly as you do. The 
whole difference is that we prove and you only assume that the process of 
iteration leads to a unique solution. Also, your so-called “energy method” 
used for the solution of your problems in elasticity and structures appears 
to me closely related to the direct methods of the calculus of variations, i.e., 
to methods which try to construct directly the minimizing function for given 
boundary values, without referring to the Euler-Lagrange differential equa- 
tion. It seems to me that after all there are many common elements in pure 
analysis and applied mathematics. 

The engineer: | shall not deny that; as a matter of fact, I have always 
felt that analysis is the backbone of applied mathematics. However, if you 
really start to apply analysis to actual cases you will see that there is a long 
way from the general idea of a method of approximation to a successful 
application of the same method. There is, for example, the question of avail- 
able time and manpower. For certain types of work, we have ingenious me- 
chanical or electrical devices such as the differential analyzer or electric 
computers. However, in most cases we have to do the computation without 
such help. Then it is not sufficient to know that the process of approxima- 
tion converges. We have to find out which method requires the least time for 
a given degree of approximation; we have to have a fair estimate of the 
improvement of accuracy by successive steps. All such practical questions 
require difficult mathematical considerations. I think we definitely need math- 
ematicians who help us to refine and, if you wish to say so, criticize and 
systematize our intuitive methods. In fact, successful applications of mathe- 
matics to engineering require the close cooperation of mathematicians and 
engineers. It is by no means a routine job to recognize the underlying common 
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mathematical relations in apparently very different fields. The mathemati- 
cian who intends to do applied mathematical research has to have a pretty 
good sense for the physical processes involved. On the other hand, the engi- 
neer has to go into the fundamentals of analysis to a considerable depth in 
order to use the mathematical tools properly. An arbitrary assembly of ma- 
chine tools does not constitute an efficient machine shop. We know there are 
powerful machine tools in your mathematical arsenal. The task before us is 
to know how to adapt and apply them. 

The mathematician: I think you’ve got something there. To carry your 
analogy further, in order to get the solution of engineering problems into pro- 
duction, you need some kind of tool designers. These are the real applied 
mathematicians. Their original backgrounds may differ; they may come from 
pure mathematics, from physics or from engineering, but their common aim 
is to “tool up” mathematics for engineering. 


A REVIEW OF THE STATISTICAL 
THEORY OF TURBULENCE* 


BY 


HUGH L. DRYDEN 
National Bureau of Standards 


1. Introduction. The irregular random motion of small fluid masses to 
which the name turbulence is given is of such complexity that there can be 
no hope of a theory which will describe in detail the velocity and pressure 
fields at every instant. Existing theories may be classified as either empirical 
or statistical. 

In the empirical theories attention is focused only on the distribution of 
mean speed and mean pressure, and assumptions are made as to the depend- 
ence of the shearing stresses required to satisfy the equations of motion of 
the mean flow. These assumptions involve one or more empirical constants. 
While the type of assumption adopted is often selected on the basis of some 
hypothesis as to the character of the fluctuations of speed and pressure, the 
theory rests on the final assumption rather than on the hypothesis as to the 
fluctuations. The various “mixing length” theories are of this type. 

In the statistical theories consideration is given to the frequency distribu- 
tion and mean values of the pressure and of the components of the velocity 
fluctuations, i.e. to the statistical properties of the fluctuations, and to the 
relation between the mean motion and these statistical properties. 

Some attempts have been made to apply the methods of statistical me- 
chanics of discrete particles. In all such attempts it is necessary to select cer- 
tain discrete elements corresponding to the particles, and to make some 
assumption as to the probability of occurrence of various values of associated 
properties or more directly the frequency distribution of the associated prop- 
erties. Difficulties are encountered at both points. The best known theory of 
this type is that of Burgers' who selected as elements in two-dimensional 
flow the points in a square network of equally spaced points and as associated 
property the value of the stream function. This theory has not as yet led to 
useful results and is not satisfactory to Burgers himself. Other attempts of 


* Received Nov. 19, 1942. 

1 Burgers, J. M., On the application of statistical mechanics to the theory of turbulent 
fluid motion, I to VII, inclusive, Verh. Kon. Akad. v. Wetensch. Amsterdam 32, 414, 643, 818 
(1929); 36, 276, 390, 487, 620 (1933). Summarized by Trubridge in Reports Phys. Soc. Lon- 
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this nature have been made by von Karman,? Noether,* Tollmien,* Gebelein,® 
Dedebant, Wehrlé and Schereschewsky,*® and Takahasi.’ 

Many of the statistical theories just mentioned do not require the turbu- 
lent fluctuations to satisfy the equations of motion nor do they require the 
fluid motion to be continuous. A statistical theory of turbulence which is 
applicable to continuous movements and which satisfies the equations of mo- 
tion was inaugurated in 1935 by Taylor’ and further developed by himself 
and by von Karman.® It is the object of this paper to give a connected ac- 
count of the present state of this particular statistical theory of turbulence. 

2. Turbulent fluctuations and the mean motion. As in other theories of 
turbulent flow, the flow is regarded as a mean motion with velocity compo- 
nents, U, V, and W, on which are superposed fluctuations of the velocity with 
components of magnitude u, v, and w at any instant. The mean values of u, v, 
and w are zero. In most cases U, V, and W are the average values at a fixed 
point over a definite period of time, although in certain problems it is more 
convenient to take averages over a selected area or within a selected volume 
at a given instant. The rules for forming mean values were stated by Rey- 
nolds’® and some further critical discussion by Burgers and others has been 
recorded in connection with a lecture by Oseen." 

When the turbulent motion is produced in a pipe by the action of a con- 
stant pressure gradient or near the surface of an object in a wind tunnel in 
which the fan is operated at a constant speed, there is considerable freedom 


2 K4rman, Th. von, Uber die Stabilitét der Laminarstrimung und die Theorie der Turbulenz, 
Proc. 1st Inter. Congr. Appl. Mech., Delft, 1924, p. 97. 

* Noether, F., Dynamische Gesichtspunkte zu einer statistischen Turbulenstheorte, Z. angew. 
Math. u. Mech. 13, 115 (1933). 

‘ Tollmien, W., Der Burgersche Phasenraum und einige Fragen der Turbulenzstatistik, 
Z. angew. Math. u. Mech. 13, 331 (1933). Brief abstract of this paper entitled, On the tur- 
bulence statistics in Burgers’ phase space, Physics, 4, 289 (1933). 

5 Gebelein, H., Turbulenz: Physikalische Statistik und Hydrodynamtk, Julius Springer, 
Berlin, 1935. 

§ Dedebant, G., Wehrlé, Ph., and Schereschewsk; 
les mouvements permanents. A pplication a la turbulence, 
(1935). Also Dedebant, G., and Wehrlé, Ph., Sur les équations aux valeurs probables d’un fluide 
turbulent, Comptes Rendus Ac. Sci. Paris 206, 1790 (1938). 

7 Takahasi, K., On the theory of turbulence, The Geophysical Magazine 10, 1 (1936). 

8 Taylor, G. I., Statistical theory of turbulence, 1-V inclusive, Proc. Roy. Soc. London 
Ser. A, 151, 421 (1935) and 156, 307 (1936). Also, The statistical theory of isotropic turbulence, 
Jour. Aeron. Sci., 4, 311 (1937). 

* K4érm4n, Th. von, On the statistical theory of turbulence, Proc. Nat. Acad. Sci. 23, 98 
(1937). Also The fundamentals of the statistical theory of turbulence, Jour. Aeron. Sci. 4, 131 
(1937). Also with Howarth, L., On the statistical theory of isotropic turbulence, Proc. Roy. Soc. 
London Ser. A, 164, 192 (1938). 

10 Reynolds, O., On the dynamical theory of incompressible viscous fluids and the determina- 
tion of the criterion, Phil. Trans. Roy. Soc. London 186, 123 (1895). 

11 Oseen, C. W., Das Turbulenzproblem, Proc. 3rd Inter. Congr. Appl. Mech., Stockholm, 


, Ph., Le maximum de probabilité dans 
Comptes Rendus Ac. Sci. Paris 200, 203 


1931, vol. 1, p. 3. 
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in selecting the time interval for which mean values are taken. So long as the 
time interval is longer than some fixed value dependent on the scale of the 
apparatus and the speed, the mean values are independent of the magnitude 
of the time interval selected and there is a clear separation between the turbu- 
lent fluctuations and the mean motion. If the mean motion itself is “slowly” 
variable, as in the case of the natural wind, difficulty arises; the separation 
becomes imperfect and arbitrary. The slowly variable mean may be taken 
over time intervals of five minutes, one day, or ten years according to the 
object of the study and the magnitude of the turbulent fluctuations varies 
accordingly. Even in flows under constant pressure gradient, there will usu- 
ally be some experimental difficulty in maintaining the conditions absolutely 
constant, and the question will naturally arise as to how the fluctuations aris- 
ing from this source may be eliminated from the “true” turbulent fluctuations. 

3. Vortex trails. For a long time every flow in which “fast” fluctuations 
of velocity occurred was regarded as a turbulent flow but experimental meas- 
urements of fluctuations show several identifiable types. The experimental 
results suggest the limitation of the term “turbulent fluctuation” to one of 
these types characterized by the random nature of the fluctuations. This 
random characteristic is in marked contrast with the regularity and periodic- 
ity noted in a second type of fluctuation associated with vortex trails. 

It is well known that when a cylinder or other object of blunt cross section 
is exposed to a fluid stream, a vortex trail appears under certain circum- 
stances, vortices breaking away with a regular periodicity. The speed fluctua- 
tions observed in the trail are periodic and in themselves do not produce 
turbulent mixing. At comparatively short distances the regular pattern trans- 
forms into an irregular turbulent motion, but the fluctuations within the trail 
itself do not have the character of the final turbulent fluctuations. 

The fluctuations of turbulence are irregular, without definite periodicity 
with time. The amplitude distribution corresponds to the Gaussian distribu- 
tion, i.e. the number of times during a long time interval that a given magni- 
tude of fluctuation is reached varies with the magnitude according to the 
“error” curve. 

If this randomness is regarded as an essential feature of the turbulent 
fluctuations, turbulence is not equivalent to any regular vortex system how- 
ever complex. The equivalent vortex picture is a large family of vortex sys- 
tems, whose statistical properties only, not individual histories, are significant. 

4. Space and time averages. The speed fluctuations u, v, and w, though 
designated the fluctuations at a point, are in reality averages throughout a 
certain volume and over a certain time as are the speed components in the 
usual hydrodynamic theory. The volume is small in comparison with the 
dimensions of interest in the flow but large enough to include many mole- 
cules. A cube of size 0.001 mm, containing at atmospheric pressure about 
2.7 X10’ molecules, satisfies this condition. The time interval is short in com- 
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parison with any time interval of interest in the mean properties of the flow 
but long in comparison with the time required for a molecule to traverse the 
mean free path. The number of collisions at atmospheric pressure is of the 
order of 5X10* per second and hence a time interval of 10-* seconds would 
suffice. 

No instruments have yet been constructed to give values averaged over 
so small a volume or so short a time interval. The best performance obtained 
to date is that of hot wire anemometers which have been developed to the 
point where average values over a cylindrical volume perhaps 0.01 mm in 
diameter and 1 mm long and over a time interval of approximately 0.5 X 10-* 
seconds can be obtained. Experimental results show that averages over these 
space and time intervals are not appreciably different from those for some- 
what larger space and time intervals and suggest that averages over smaller 
intervals would not be appreciably different. The results also suggest that 
measuring equipment that does not approach these space and time intervals 
gives results which largely reflect the properties of the measuring instrument 
rather than the properties of the turbulent fluctuations. In other words the 
measurement is that of a variable mean velocity over space and time intervals 
fixed by the characteristics of the instrument, rather than measurements of 
the turbulent fluctuations. If the frequency spectrum of the turbulent fluctua- 
tions is known, the effect of the instrument characteristics can be estimated, 
as discussed in section 19. 

5. Pulsations. Reference has previously been made to the difficulty in 
certain cases of making a clear separation between the mean motion and the 
turbulent fluctuations, because of the difficulty of defining a time interval 
long enough to include many fluctuations but small enough so that the mean 
varies only slowly. The difficulty is often increased by the presence of a fairly 
rapid variation of the mean speed over large areas, perhaps the entire cross 
section of the fluid stream, to which the name pulsation may be given. Sucha 
fluctuation is recognizable by the fact that there is a regularity in the space 
distribution of the fluctuations such that definite phase relations exist. Pulsa- 
tions have been observed in laminar flow in boundary layers. An essential 
characteristic of the turbulent fluctuations is an irregularity and randomness 
in the space distribution as well as in the time distribution. 

It is often possible to eliminate the effect of pulsations on the measure- 
ments by a low frequency cut-off in the equipment for measuring u, v, and w. 
The choice of the cut-off frequency is equivalent to a selection of the time 
interval over which averages are taken to obtain the mean speed and by this 
device the pulsations are regarded as variations of the mean speed. 

6. Continuity of the turbulent motion. It is well known that the structure 
of a fluid is in the final analysis discontinuous, the fluid consisting of individ- 
ual molecules. Nevertheless the usual hydrodynamic theory regards the fluid 
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as a continuum. Such an assumption can be justified when the dimensions 
of the flow system are very large compared to the mean free path of the 
molecules. The velocity of the fluid at any point is then defined as the vector 
average of the velocities of the molecules in a small volume surrounding the 
point, the value obtained being independent of the magnitude and shape of 
the volume within certain limits. 

Some investigators” have concluded that the phenomena of turbulence 
require the assumption of discontinuity in the instantaneous components. 
The Taylor-von Karman statistical theory retains the assumption that the 
fluctuations are continuous functions of space and time as in Reynolds’ theory. 

The applicability of this assumption is a matter for experimental determi- 
nation. If experimentally a volume and time interval can be selected which 
may be regarded as large in comparison with molecular distances and periods 
but small as compared to the volumes and time intervals of interest in the 
turbulent fluctuations, the fluctuations may be safely regarded as continuous. 
As described in section 4, the experimental data perhaps do not prove but do 
definitely suggest that such a choice is possible and to that extent the as- 
sumption of continuity is experimentally justified. 

7. The Reynolds stresses. If in the Navier-Stokes equations of motion the 
components of the velocity are written as U+-u, V-+v, W+w, thus regarding 
the motion as a mean motion U, V, W, with fluctuations u, v, w superposed, 
and mean values taken in accordance with the rules mentioned in section 2, 
a new set of equations is obtained which differs from the first only in the pres- 
ence of additional terms added to the mean values of the stresses due to vis- 
cosity. These additional terms are called the Reynolds stresses or eddy 
stresses. The eddy normal stress components are —pu?, —pv?, —pw* and 
the eddy shearing stress components are —puv, —pvw, —puw. Each stress 
component is thus equal to the rate of transfer of momentum across the cor- 
responding surface by the fluctuations. 

In the light of kinetic theory the eddy stresses closely parallel in origin 
the viscous stresses. It has been explained how u, v, and w are themselves the 
mean speeds of many molecules. The effect of the molecular motions appears 
in the smoothed equations of the continuum as a stress, the components of 
which are equal to the rate of transfer of momentum by the molecules across 
the corresponding surfaces. 

8. Correlation. If the fluctuations were perfectly random, the eddy shear- 
ing stress components —puv, —pvw, —puw would be zero. The existence of 
eddy shearing stresses is dependent on the existence of a correlation between 
the several components of the velocity fluctuation at any given point. The 
coefficient of correlation between u and 7 is defined as 


7 


2 Kampé de Fériet, J., Some recent researches on turbulence, Proc. Fifth Inter. Congr. 
Appl. Mech., Cambridge, Mass.,1938, p. 352. 
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uD 
R.. = — =" (8.1) 
V urrv/v? 


The mean values 22, \/v?, and \/wz? are often called the components 
of the intensity of the fluctuations. 

The eddy shearing stress may be written in terms of the correlation coeff- 
cient as 

— puv = — pRuvV/uvv? (8.2) 
and similarly for the other components. 

In addition to the correlation between the components of the velocity 
fluctuations at a given point, the Taylor-von Karman theory makes much 
use of correlations between the components of the velocity fluctuations at 
neighboring points. Denote the components of the fluctuations at one point 
by m1, 1, %:, and at another point by we, ve, we. The coefficient of correlation 


between 2; and 72 is defined as 





Rey * ae (8.3) 


9 


and similarly for any other pair. These correlation coefficients form useful 
tools to describe the statistical properties of the fluctuations with respect to 
their spatial distribution and phase relationships. 

9. Scale of turbulence. The earliest attempt to describe the spatial char- 
acteristics of turbulence was the introduction of the mixing length concept, 
the mixing length being analogous to the mean free path of the kinetic theory 
of gases. Logical difficulties arise because there are no discrete fluid particles 
in the turbulent flow which retain their identity. A method of avoiding these 
difficulties was suggested by Taylor many years ago. He showed that the 
diffusion of particles starting from a point depends on the correlation R, be- 
tween the velocity of a fluid particle at any instant and that of the same par- 
ticle after a time interval ¢. If the functional relationship between R; and ¢ is 
of such a character that R; falls to zero at some interval T and remains so for 


greater intervals, it is possible to define a length /; by the relation: 


T “ 
h=y, # Ridt = Ff R,dt (9.1) 
0 0 


in which v is the component of the velocity fluctuations transverse to the 
mean flow and in the direction in which the diffusion is studied. 


3 Taylor, G. I., Diffusion by continuous movements, Proc. London Math. Soc. Ser. A, 20, 
196 (1921). 
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This method of assigning a scale to turbulence is of value in the study of 
diffusion as described in section 22. It is based on the Lagrangian manner of 
describing the flow by following the paths of fluid particles. It is more com- 
mon to use the Eulerian description by considering the stream lines existing 
in space at any instant. Taylor later® suggested a method of describing the 
scale in the Eulerian system based on the variation of the correlation coeffi- 
cient Ry between the values of the component wu at two points, separated by 
the distance y in the direction of the y coordinate, as y is varied. The curve 
of R, against y represents the statistical distribution of u along the y axis at 
any instant. If R, falls to zero and remains zero, a length L may be defined 


by the relation: 


L= f R,dy. (9.2) 
0 


The length Z is considered a possible definition of the average size of the 
eddies present and has been found to be a most useful measure of the scale 
of the turbulence, especially for the case of isotropic turbulence. Correspond- 
ingly, a length L, may be defined by the relation: 


t= f R.dx (9.3) 
0 


where R, is the correlation between the values of the component u at two 
points separated by distance x in the direction of the x coordinate. 

10. Isotropic turbulence. The simplest type of turbulence for theoretical 
or experimental investigation is that in which the intensity components in 
all directions are equal. More accurately, isotropic turbulence is defined by 
the condition that the mean value of any function of the velocity components 
and their derivatives at a given point is independent of rotation and reflection 
of the axes of reference. Changes in direction and magnitude of the fluctua- 
tions at a given point are wholly random and there is no correlation between 
the components of the fluctuations in different directions. Thus 42=0?=w* 
and uv=vw = uw =0. 

There is a strong tendency toward isotropy in all turbuient motions. The 
turbulence at the center of a pipe in which the flow is eddying or in the natu- 
ral wind at a sufficient height above the ground is approximately isotropic. 
A grid of round wires placed in a uniform fluid stream sets up a more or less 
regular eddy system of non-isotropic character which very quickly transforms 
into a field of uniformly distributed isotropic turbulence. 

The assumption of isotropy introduces many simplifications in the statis- 
tical representation of turbulence. The two quantities, intensity and scale, 
appear to give a description of the statistical properties of the turbulent field 
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which is sufficient for most purposes. Turbulent fields of this type can readily 
be produced experimentally and studied. The intensity may be varied from 
less than 0.1 to about 5.0 percent of the mean speed and the scale independ- 
ently from a few mm to 25 mm." 

11. Decay of isotropic turbulence. The kinetic energy of the turbulent 
fluctuations per unit volume is equal to 3o(u?+-v?+w*) which for isotropic tur- 
bulence becomes (3/2)pu?. The rate of decay is therefore —(3/2)pd(u?) /dt. If 
the isotropic turbulence is superposed on a stream of uniform speed U, we 
may write dt=dx/U and hence the rate of decay with respect to distance x 
as —(3/2)pUd(u?) /dx. 

In a fully developed turbulent flow the Reynolds stresses are proportional 
to the squares of the turbulent fluctuations. The work done against these 
stresses, which in the absence of external forces must come from the kinetic 
energy of the system, is proportional to pu’*/L where uw’ is written for £7 
and L is a linear dimension defining the scale of the system, which may be 
taken as the LZ defined by (9.2). Equating the two expressions for the dis- 
sipation and designating the constant of proportionality as 3 A, we find: 


— (3/2)pUd(u'?)/dx = 3Apu’?/L (11.1) 
or 
Ld(U/u')/dx = A. (11.2) 
Integrating: 
U/u’ — U/us = af dx/L (11.3) 


where U/ug is the value of U/u’ at x=x». This equation has been found to 
give a very good representation of the experimental data. The essential fea- 
tures of the derivation were given by Taylor. To evaluate the integral, L must 
be known as a function of x. Taylor’s first proposal was to assume that L is 
independent of x and proportional to the mesh M of the grid giving rise to 
the turbulence. If Z is constant, 


U/ul — U/ug = A(x — x0)/L (11.4) 


giving a linear variation of U/u’ with x. Assuming L/M=k, Taylor found 
values of A/k for data from various sources varying between 1.03 and 1.32. 


4 Dryden, H. L., Schubauer, G. B., Mock, W. C., Jr., and Skramstad, H. K., Measure- 
ments of intensity and scale of wind-tunnel turbulence and their relation to the critical Rey- 
nolds number of spheres, Tech. Rept. Nat. Adv. Comm. Aeron. No. 581 (1937). 
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When measured values of L became available it was found that L in- 
creased as x increased, the results being represented empirically within the ac- 
curacy of the measurements by the relation ZL =Lo+c(x—xo), whence 


) 


wm 


U/u' — U/us = (A/c) log. [1 + c(x — x0)/Lo]. (11. 


Taylor found values of A for data from various sources varying between 
0.43 and 0.19. 

Further study suggests another relation for the variation of Z with x. 
A discussion of the general theory will be deferred until section 17 and the 
question discussed on purely dimensional considerations. If one assumes that 
du’ /di, the rate of change of intensity, and dL /dt, the rate of change of scale, 
are determined solely by the values of Z and w’, i.e. that viscosity and up- 
stream conditions have no influence, it follows from dimensional reasoning 


that 


B (11.6) 


Ld(1/u’)/dt = A and (1/u’)dL/adt 


or 


B (11.7) 


Ld(U/u')/dx = A and (U/u’)dL/dx 


where A and B are numerical constants. The first equation of each pair is 
the same as equation (11.2); the second is a new relation. 
Integration of equations (11.6) and (11.7) leads to the relations: 








ug r (A + B)ud (x — x9) 414+) 
e138. fee (11.8) 
v 6k LU t 
and . 
L q (A + B)ug (x — x9) ]2/4+®) 
=jit 1" en (11.9) 
Lo L LoU a 


where ug and Lo are the values at x =0. 

[f it is desired to introduce a reference dimension pertaining to the dimen- 
sions of the grid producing the disturbance, this may be done, but according 
to equations (11.8) and (11.9) any dimension may be used and the decay 
does not depend on its value. The mesh distance M is often used but certain 
results reported by von Kaérman" show that if 1//d is not too small, the use 


6 Taylor, G. I., Some recent developments in the study of turbulence, Proc. Fifth Inter. 
Congr. Appl. Mech., Cambridge, Mass., 294 (1938). See later detailed report of measurements 
in Hall, A. A., Measurements of the intensity and scale of turbulence, Rept. and Memo. No. 
1842, Aeronautical Research Committee, Great Britain (1938). 

16 K4rm4n, Th. von, Some remarks on the statistical theory of turbulence, Proc. Fifth Inter. 
Congr. Appl. Mech., Cambridge, Mass., 1938, p. 347. The grid dimensions are not given in 
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Fic. 1. The turbulent fluctuation u’ behind a grid of wires of diameter d 
as a function of distance x from the grid. 
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Fic. 2, The scale ZL behind a grid of wires of diameter d as a 
function of distance x from the grid. 
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of the wire diameter d as the reference dimension leads to a single curve for 
all grids irrespective of the mesh-diameter ratio. 

The available data are plotted in Figs. 1 and 2 from references in foot- 
notes 14, 15, and 16. The solid curves are respectively 


(U/u’)? = 400[(1 + 0.04(«/d — 80) } (11.10) 


to 


and 
(L/d)? = 0.264[(1 + 0.04(x/d — 80) | (11.11) 
which are in the form of equations (11.8) and (11.9) with the constants 
A =B=0.2056. These curves are frankly selected to fit the National Bureau 
of Standards data. 

If one considers the complete system of screen and turbulent field, dimen- 
sional considerations suggest that for geometrically similar screens whose 
scale is fixed by some characteristic dimension, such as the mesh length M, 
the ratios u’/U and L/M would be a function of x/M, of the Reynolds Num- 
ber UM/v and of the turbulence of the free stream u/ /U, in which the screen 
is placed. If the screens are not geometrically similar but are made up of 
cylindrical rods of diameter d, the intensity and scale also depend on d/M 
and on the roughness of the screen. The effects of these parameters have not 
been fully investigated, and doubtless a part of the discrepancy between the 
available results is to be ascribed to the influence of these factors. 

For example, the screens used at the National Bureau of Standards were 
either woven wire screens or wooden screens with fairly rough surfaces with 
the members interlacing in the wire screens and intersecting in the wooden 
screens. The ratio d/M varied from 0.186 to 0.201. The screens used by Hall 
were arranged in two planes, i.e., horizontal rods in one plane, vertical rods 
just touching the horizontal rods but in another plane. The ratio d/M was 
0.184 to 0.188. Von Karman has studied the effect of varying d/ M from 0.086 
to 0.462 and has used screens both of the woven type (results published by 
von Karman, loc. cit.) and of the biplane type (results not published). A 
study of these data suggests that the difference between the results for woven 
screens and biplane screens is unimportant and that if results are plotted in 
terms of x/d rather than x/M the effect of d/M is small for values of d/M 
near 0.2. No data are available on the effect of roughness. 

Few data are available on the effect of free stream turbulence. Hall ob- 
tained an increase of about 10 to 20 percent in wu’ for a 1-inch screen at the 








the paper, but Professor von K4rm4n has kindly supplied them as follows: 


Grid Mesh Dista nce, M Wire Diameter, d M/d 
inches inches 
1 4.96 0.230 2.16 
2 5.00 - 105 4.75 
3 5.07 .084 6.03 
4 4.99 043 11.6 
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same value of x/M by increasing the free stream turbulence from 0.2 percent 
to 1.3 percent. We have had the opportunity of making some measurements 
behind the same 1-inch screen used in the measurements described in NACA 
Technical Report No. 581 in an airstream for which the free stream turbu- 
lence is 0.03 percent as compared with 0.85 percent for the older measure- 
ments. The results are shown in Fig. 3 as compared with Hall’s measurements. 
It is obvious that the turbulence of the free stream is one of the controlling 
factors, but not the only one. 
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Fic. 3. Effect of free stream turbulence on the turbulence behind a 1-inch screen. 


The study of the turbulent field behind screens as affected by numerous 
parameters is of interest from the standpoint of a study of screens. However, 
the turbulent field may be regarded from another point of view, i.e. in relation 
solely to the theory of isotropic turbulence. If the turbulence is truly iso- 
tropic, and if its characteristics can be adequately described by the two 
quantities, intensity and scale, its behavior can depend only on the values 
of intensity and scale at some given point. The details of construction of the 
source screen and its distance upstream are of no importance. Even the in- 
fluence of the turbulence of the free stream should be absorbed in the given 
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values of u’ and L at some one point. The decay of isotropic turbulence is 
considered from this point of view in section 17. 

12. Effect of contraction. The behavior of turbulence in a contracting 
stream is of interest in connection with the flow in the entrance cone of a 
wind tunnel. Prandtl!’ suggested that the longitudinal components of the 
fluctuations were reduced in the ratio of 1 to / where / is the ratio of the en- 
trance area to the exit area of the cone. This result was derived on the as- 
sumption that the gain in energy is the same for all filaments traversing the 
cone. The same result was obtained from the Helmholtz vortex theorem, 
which was also used to show that the lateral components were increased in 
the ratio »//. Since the mean speed increases proportional to /, the values of 
u’/U and v’/U are reduced according to this theory in the ratios 1//? and 1/+/1 
respectively. This computation neglects the decay of the turbulence because 
of viscosity. 

Taylor'* computed the effect of a contraction on certain mathematically 
defined forms of disturbance. Two objections may be offered to this treat- 
ment. First, as in Prandtl’s treatment, the decay of the turbulence is neg- 
lected. Second, the computation is made on a regular disturbance which is 
assumed to retain its regularity. When the rapid development of an isotropic 
turbulent field from a Karm4n vortex trail is considered, it is hard to believe 
that a regular vortex pattern could retain its character throughout the length 
of a wind tunnel entrance cone unless the scale was very large indeed. 

If it is assumed that the istropic turbulent field is unaffected by changes 
in the mean speed, the decrease in u’ may be computed from the decay during 
the time required for the fluid to traverse the cone. This time interval is 
Jidx/U. If A is the area of the cross section at any value of x, UA =U Ao 
where U» and Ao are the values at x =xo, and hence the time interval is 
JRA dx/UvA. 

There are as yet no suitable experimental data for checking any theory. 
In the measurements quoted by Taylor all the data were obtained sufficiently 
close to a grid to lie within the non-isotropic turbulence of the vortex trails 
from the individual wires. 

13. The correlation tensor function. Von Karman’ introduced the cor- 
relation tensor function in the statistical theory of turbulence as a generaliza- 
tion of the particular correlation coefficients discussed by Taylor. The cor- 
relation coefficients between any component of the speed fluctuation at a 
given point and any component of the speed fluctuation at another point 





17 Prandtl, L., Herstellung einwandfreier Luftstréme (Windkandle), Handbuch der Experi- 
mentalphysik, F. A. Barth, Leipzig, 1932, Vol. 4, Part 2, p. 73. 

18 Taylor, G. I., Turbulence in a contracting stream, Z. angew. Math. u. Mech. 15, 91 
(1935). 

19 K4rm4n, Th. von, and Howarth, L., On the statistical theory of isotropic turbulence, 
Proc. Roy. Soc. London, Ser. A, 164, 192 (1938). 
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form a tensor. If one point is held fixed and the other varied, the tensor varies 
as a function of the coordinates of the variable point with respect to the fixed 
point. We may speak of this function as the correlation tensor function. 
In isotropic turbulence the correlation tensor has spherical symmetry and 
the several components are functions only of the distance r between the two 
points, and of the time ¢. Denote by m, 1, #1 and ue, v2, We the components 
of the velocity fluctuations at the two points P; and P2 having coordinates 


(x1, 0, 0) and (x2, 0, 0) respectively. Suppose that uj, vj, wj, which by isotropy 

are equal, are independent of position and equal to 43. Then 43 =13 =u =u’. 
The correlation coefficients 1,22/u? and w,w2/u? will be identical because 

of isotropy and will be some particular function of the distance r between 


Z 








Y 


Fic. 4. 


P; and P, and of the time ¢, say g(r, #). The correlation coefficient 1,u2/u? 
will also be a function of r and ¢t, say f(r, t). The correlation coefficients 
UyV2/ U2, UjWe/ U2, V,U2/U2, V{We/ U2, W{U2/Uu? and w,v2/u? can be shown to be zero. 
Thus if the Y and Z axes are rotated about the X axis through 180°, the ab- 
solute values of all components are unchanged but the signs of the v and w 
components are reversed. Denoting values referred to the new axes by capital 
letters, Uj=m, U2=u2, Vi=—1, Ve= —v2, Wi=wi, We= —we, so that, for 
example, U; V2= — 1%. But by isétropy, the value of any function of the 
components is unchanged by rotation of the axes, and therefore U; Vo = uy02. 
To satisfy both relations u,v, must equal zero. Similarly for the other terms 
containing or u». By reflection in the XZ plane v,w2 and w,¥2 may likewise 


be shown to be zero. 
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The correlation coefficient for components of the fluctuations in any arbi- 
trary directions at any two points may be expressed in terms of the functions 
f(r, t) and g(r, t) and the geometrical 

parameters. Consider any two points 

> > P and Q and components of the fluc- 
eee . nee tuations p in the direction PP’ at P 
and q in the direction QQ’ at Q. (Fig. 

4). Denote by QQ” the orthogonal 


u, projection of QQ’ on the plane PP’Q; 
by a, 8B, and y the angles P’PQ, 
7 — PQQ”, and QQ’Q”’; and by fi, po, 


ps and qi, G2, gs the components of the 
fluctuations at P and Q in the direc- 


ths F(a,t) U, 











ae (n,t) See 
ii tion PQ, in the direction normal to 
PQ and Q’Q”, and in the direction 
Q’’Q’. Then 
p = picosa+ pesina 
| Uy ee es 
, g = qi cos B sin y + gz Sin B sin ¥ 
Fic. 5. The principal double correlations + qs COS ¥. (13.1) 
in isotropic turbulence. 
Hence 
Pq = pigi cos a cos B sin y + pog2 sin B sin a sin (43.2) 


the other terms vanishing as proved in the preceding paragraph. In terms 
of f(r, ¢) and g(r, t) 


Pg, ‘“e = f(r, t) cos a cos 8 + g(r, t) sin @ sin B} sin ¥. (13.3) 


The correlations denoted by f(r, ¢) and g(r, ¢) are indicated in Fig. 5. 
If now any two points with coordinates (x1, 1, 21) and (x2, ye, 22) and speed 
fluctuations with components 1, v1, w; and uo, v2, We are considered, the nine 





quantities UjU2, UyV2, UyWe, VyUy, V{V2, VjWe, WU, WyV2, and w ,We2 are the compo- 
nents of a second rank tensor. Each one may be evaluated by equation (13.2) 


Pi 


in terms of f(r, ¢) and g(r, ¢) with the result in tensor notation 


f(r, t) — g(r, t) 
R ee dy Sl ik t)I (13.4) 


9 


72 
where r is the vector having components X =x2—x1, Y=ye—-y1, Z=22—21 
100 


ris |r| and Jis the unit tensor |0 1 0 
(001) 
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The velocity fluctuations satisfy the equation of continuity. Hence 


Ou2 Ove OW2 








= 0. (13.5) 





OX2 OV2 O22 
Multiplying by u:/u? which is independent of x2, ys, 2 and introducing the 
correlation coefficients Ry,u,, etc. and the components X, Y, Z of r: 


IRu jus 4 ORu, vs Pip, 0 (13 6) 
ox oY 0Z ii . 








From equation (13.4) 


Boe bt “=< 
= ano A* + g; Raw = - : as: Ru ws as f 
r2 


r? r? 


XZ 








whence, remembering that X?+ ¥°+2Z?=r?, dr/8@X =X/r, Of/0X =(df/dr) 
(dr/0X) =(X/r)(Of/dr), etc., equation (13.6) becomes 


X[2(f — g) + r(af/dr)] = 0. (13.7) 
The continuity equation must be true for any value of X. Hence 
2f(r, 4) — 2g(r, 1) = — rdf(r, t)/or. (13.8) 


The correlation tensor can thus be expressed in terms of a single scalar 
function, either f(r, 2) or g(r, t). The function g(r, ¢) is the correlation coeff- 
cient previously denoted by Ry. The scale L=f-R, dy=f-g dr. The integral 
ScR: dx =f dr is termed the longitudinal scale L, to distinguish it from 
the lateral scale L. Obviously from equation (13.8) 


L-L.= rf r(0f/dr)dr = if x(OR./dx)dx. (13.9) 
0 0 
Since f and g are even functions of r, 
j= 1+ for?/2+--- (13.10) 
g=1+g%r2/2+.---. (13.11) 


From equation (13.8), 2f6 =g9, whence for small values of r, 
[1 + (6/2)r?]2 + [(fo — g6/2]rr 

(1 + for)I — (1/2)ffer. 

We require later the second derivatives of R at r=0, ie. X= Y=Z=0, as 
follows: 

0*R, ue? FPR... P* Rew 

eee we wes >= f% (13.13) 


0X? oY? 0Z? 


R 


(13.12) 
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O* Ruy,  O?Ruyus as P 
= ——— and similar terms obtained-by cyclic exchange = 2f¢ (13.14) 





ay? QZ? 
8 Ruive ° ° ° ° . ~ 
rong and similar terms obtained by cyclic exchange = — (1/2) /f¢. (13.15) 
AC 
O° Rusvs 





All others, e.g. etc. are zero. (13.16) 


0Xd0Z 

Von Kérman points out that the correlation tensor is of the same form 
as the stress tensor for a continuous medium when there is spherical sym- 
metry. In the analogy f(r) corresponds to the principal radial stress at any 
point, g(r) to the principal transverse stress, and the several R’s to the stress 
components over planes normal to the coordinate axes. The relation between 
f and g given by the continuity equation corresponds to the condition for 
equilibrium of the stresses. 

Equation (13.8) has been experimentally checked at the National Physi- 
cal Laboratory. 

14. Correlation between derivatives of the velocity fluctuations. In fur- 
ther developments it will be necessary to know the mean values of the prod- 
ucts of the derivatives of the components of the fluctuations at a given point, 
for example (8u;/9x;) (80,/8y,). These mean values may readily be computed 
from the correlation tensor. Thus: 





8(u302) a O(Ru,: ) — ARe. ) 
bith bes Pagp petite. en pa =. (14.1) 
Oxy Ox ax 





Since v2 is not a function of x1, this may be written 











(0u;/0x1)02 = — u?OR,,»,/OX. (14.2) 
Differentiating now with respect to y2 
Ou, OX» 0 fou, — 9 /ORa», —. Ry, 0, 
—  — = —| — 2} = — 4? —| — ] = — uv? ——_-- (14.3) 
0X1 O42 OV2\Ox, Oe ox s aX0Y 
Now letting P; and P: coincide, 
Ou dv 07 Ruv\ : 
——_—_—_— = u? s “ x => } = 0. (14.4) 
dx dy oXd0Y 


The limiting value of the second derivative has previously been computed 
(equation (13.15)), whence 


F wn (14.5) 
x Oy 2 


By similar reasoning it may be shown that 
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uA? Jai? fw 
(“) a (=) -( ) = - aif (14.6) 
Ox oy Oz 


(=) - (2) -(&) 
ax) . \as/ \ax 








(=) (“*) 
oy Oz 


(14.7) 


(aw? vm 
—]} = — 2u*fo 
oy 


du Ow 


II 


and 
dv Ou dw dv “ 
——=— -—= — = — ft. (14.8) 
Ox dy Oy Oz Oz Ox 2 
The method can be extended to derivatives of higher order. 
15. Triple correlations. Von Kar- 
man designates the mean values of vi h(n,t) Ue 
the product of three components 





— 
of the velocity fluctuations, two of ” 
which are taken at one arbitrary 
point and the third at a second arbi- u's K(a,t) - 

- 





trary point, as triple correlations. 
They arise when correlation coeff- 
cients are introduced into the equa- 
tions of motion. He shows that the | “Vy 
triple correlations are components of > 
a tensor of third rank designated T “2 
which is a function of X, Y, Z and the 
time. He proves that in isotropic tur- rr. 
bulence this tensor can be expressed oi 7 ; , 
. . Fic. 6. The principal triple correlations 
in terms of three functions h(r, ¢), in isotropic turbulence. 
k(r, t), and g(r, t) corresponding to 
the correlations shown in Fig. 6, and that the development of these func- 
tions in powers of r begins with the r° term. The equation of continuity 
permits the expression of k and gq in terms of h by the relations: 
k= —2h (15.1) 
q = — h — (r/2)(dh/dr). (15 .2) 
Thus the tensor T can be expressed in terms of a single scalar function h(r, 2). 
16. Propagation of the correlation with time. The fluctuations are assumed 
to satisfy the equations of motion, namely, 


q lat) «| 








Ou, Ou Ou, Ou, 
— + + +w 


1 
Ox Oy) 02 


1 Op O7u O7u 07 uy 
a ey ee ee + a ~) (16.1) 
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and the two equations obtained by cyclic permutation. 
Multiplying this equation by we, introducing X, Y, and Z, and taking 


mean values: 














Ou,  O(uiue) O(uy2U2) (4 Wyte) 
ue — — ———- —- ——__ —- ———_ 
ot ox Oo} OZ 
1 Op 02 usu» 02 \Us O07 U2 
= — — u— + »( —— + ——_ + ——}._ (16.2) 
p Ox, 0x? oY? oZ? 


By an analogous procedure, it may be shown that: 





Ouse Ousuy OUgVoy OU2Wot, 
“y— — —— —- ———_ - —— 
Ot OX oY OZ 
— a ‘ — (16.3) 
1 Op O07; Ue O7u\Ue O07 uu 
ae gre +>( Poe o — a se 
p OX OA* oY? 0Z? 


Von Karman shows that the pressure terms vanish. Adding the two equations 
and introducing the correlation coefficients, we find 


a <= 0 a as a oO Ss enh ee 
— (#*Ru,u,) — (4?)®/? — (ugue + ugui) — (u*)®!? — (uj01%2 + UVM) 
ot ox oy 
is 0 a a ae 
9\2/9 ‘ 
ee ee (U1 WU. + U2Wett;) (16.4) 
d, 


—_[ 0?u1u2 O07 uu O07 ule 
ht ——- 44+ — I. 
0X? oY? 0Z? 
This equation may be expressed in terms of the functions f, g, k, g and A. 
Then by using the relations (15.1) and (15.2) between these functions ob- 


tained from the equation of continuity, a partial differential equation between 


f and k is obtained, namely 


(fu?) — (doh th ee es | 4 of , 
— + 2a (= + - ) = 20a ( +— - ). (16.5) 


Ot Or r Or? r or 


This is the equation for the change of the function f with time, but it can- 
not be solved without some knowledge of the function h. 

17. Self-preserving correlation functions. Let us suppose that the func- 
tions f(r, ¢) and h(r, ¢) preserve the same form as ¢ increases, only the scale 
varying. Such functions will be termed “self-preserving.” If ZL is some measure 
of the scale of the correlation curve, f and h will be functions of r/Z only, 
where L is a function of t. The length Z may be any measure of the scale such 
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as the radius of curvature of the correlation curve at r=0 or any other desig- 
nated point, the value of r for a given value of the correlation coefficient, or 
the quantity obtained by integration of the correlation coefficient from r=0 
to infinity which has previously been termed the scale of the turbulence. 
Introducing the new variable y=r/L and placing (u?)”?=u’ in equation 
(16.5), we obtain 


fL dw? 1 dL af ah hy 2/af 44 
aloo 2 42(— 442) -—(S 425) ary 
u® dt u’ dt ody Oy y N \oy? y oy 





where N is the Reynolds Number of the turbulence u’L/v. Since the coeffi- 
cient of the third term is a numerical constant, the functions f and h will be 
functions of ¥ and ¢ alone only if the coefficients of the other terms are also 
numerical constants. This requires that 





L du’? d(1/u’ 
bd (17.2) 
u's dt at 
1 dL 
inate sien te (17.3) 
u’ dt 
u'L 
oe 7 (17.4) 
Vv 


where A, B, and No are independent of u’, L, and ¢. It is readily shown that 
these relations are consistent only if A =B and that the solutions are 





1 1 2A 
slams sai saaed a t (17.5) 
u’? ug? N ov 

i -f,< tite (17.6) 


where u¢ and Lo are the values for ¢=0 and ug Lo/v = No. 

These equations are in the form of equations (11.8) and (11.9) with A=B 
and agree well with the formulation of the experimental data represented by 
equations (11.10) and (11.11) with 2A Nov/L? =2Auj/Nov=2A V/ui/Ly = 0.04, 
corresponding to ¢=0 at a distance of 80 wire diameters from the grid. The 
constant A is equal to 0.2056 when L is defined as S-Riay. 

For self-preserving turbulence equation (16.5) becomes 


—Af—Ay(af/dap)+2(dh/apy+4h/p) = (2/No) [(02f/ay?+(4/y)(af/ay)]. (17.7) 


This equation determines the shape of the correlation curve. Von Kérman’® 
discusses the shape when the function h is neglected. The shape depends on 
the Reynolds Number Np of the turbulence. The shape also depends on the 
constant A but closer examination shows that A is always associated with L 
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and is dependent on the method of defining L. If y is set equal to rA/L in- 
stead of r/L, and the length L/A is used instead of L in the definition of the 
Reynolds Number of the turbulence, the A disappears from equation (17.7). 
Whether the values of L defined by /gdr will yield the same values of A for all 
shapes of correlation curves described by (17.7) cannot be definitely an- 
‘swered. 

Approximate solutions of (17.7) are not easy since it turns out that f 
varies with No in such a manner that, for small values of at least, the term 
on the right-hand side is of the order of unity. 
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Fic. 7. Hall’s measurements of turbulence behind screens. 


According to this suggested theory, the shape is self-preserving and the 
Reynolds Number remains constant during the decay of a given turbulent 
field. The scale approaches very large values as the’intensity approaches very 
small values. The length A (which is discussed in section 18) is proportional 
to L. For different values of the Reynolds Number of the turbulence the con- 
stant of proportionality varies inversely as the square root of the Reynolds 
Number. Likewise the shape of the correlation curve varies with the Reynolds 
number of the turbulence. 

Equation (17.6) shows the same functional relation between the scale and 
the time as given by Prandtl at the Turbulence Symposium as a result of his 
analysis of photographs of the decay of isotropic turbulence. 
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Von Karman also discusses the case in which the assumption is made 
that the self-preserving feature applies only to large values of y and the 
Reynolds Number N> is sufficiently large that the right-hand term of (17.7) 
can be neglected. In this case (17.2) and (17.3) are obtained without (17.4) 
and the solution is identical with that given by (11.8) and (11.9) of sec- 
tion 11. The theoretical equations (17.5) and (17.6) do not involve either U 



































16 1 ) 
| | | ° 
4). “SCREEN No2| | 7 
+ SCREEN No.3 | | 
12 | 0 SCREEN No.4 | 
| 
10 }— —_ 
u) 8 | 
6 — - —— 
| . i 
a wr 
pet | 
2 ; Lr a eS — oe 
ot | | 
pe | it | | 
< 0 2 4 6 8 io 2 14 
ut 
= 


Fic. 8. Von Karmén’s measurements of turbulence behind screens. 


or M explicitly. However, for comparison with experimental data, t 2y may 


be written as follows: 


uo \? 2Augt uo (xX — Xo) ' 

( )=14 -=1+2A — —— (17.8) 
u’ Lo l Lo 
L\? 2Augt ug (x — Xo) 

( )=1+ —— = +24 — : (17.9) 
Lo Lo l Lo 


Both ué and Lo should be known, but unfortunately Z» was not measured 
in all of the experiments. 

Figures 7, 8, and 9 show the results ot Hall, of von Karman, and of the 
author and his associates (designated NBS) plotted in a manner to facilitate 
comparison with equation (17.8). 

The reference position x» has been taken as 40 times the mesh length ex- 
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cept for von KérmAn's results for which xo was taken as 212.5 times the rod 
diameter (equivalent to xo/ M =40 for d/ M =0.188). In the absence of definite 
information as to Lo, Lo/M was assumed equal to 0.29 except for von Karman’s 
results for which Lo/d was assumed to be 1.54 (equivalent to Lo/M=0.29 for 
d/M=0.188). The value of u@ was determined by interpolation from the ob- 
servations of each experimenter near x/M=40, giving the following results. 
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Fic. 9. NBS measurements of turbulence behind screens. 
Rod ‘ 
P A esh - d Spee F 
Experimenter _— Diameter - pur ug /l Remarks 
Inches ft/sec 
Inches : 
Hall 1.0 0.188 20 0.0146 
0.5 .092 20 .0144 
0.5 .092 40 .0152 
0.5 .092 80 .0174 
von K4rm4n 0.5 0.105 38 and 54 .0201 Screen 2 
0.5 .084 38 and 75 .0201 Screen 3 
0.5 .043 38 and 75 .0299 Screen 4 
NBS 0.25 0.050 20-70 0.0250 NACA Tech. Rept. 581 
0.5 .096 20-70 .0221 NACA Tech. Rept. 581 
1.0 .196 20-70 .0224 NACA Tech. Rept. 581 
1.0 0.196 30 0.0188 Recent tests 


1.0 .196 70 .0173 Recent tests 


1943] A REVIEW OF THE STATISTICAL THEORY OF TURBULENCE 31 


The values of ug /U range from 0.0144 to 0.0299; presumably the differ- 
ences are due mainly to the factors discussed in section 11, although system- 
atic errors may be partly responsible. 

In each figure, equation (17.8) with constant A equal to 0.29 is plotted 
as a straight line. Most of the points would be better fitted by a curve of 
increasing slope with increasing time. It thus appears that equation (11.9) 
with (A+B)/A having some value between 1 and 2 fits the experimental data 
better than (17.8). 

However, the data are not at all consistent. The departures are largest 
for the smaller values of u’/U. In Hall’s experiments, the results on the 
+-inch screen show little systematic departure at 40 and 80 ft/sec, whereas 
those on the same screen at 20 ft/sec and on the 1-inch screen begin to rise 
above the line at ug t/L)o=0.5. Von Karman’s data on screen 2 at 38 ft/sec 
lie near the line; those on the same screen at 54 ft/sec and on screens 3 and 4 
at 38 and 75 ft/sec begin to rise above the line at uj t/Lo=4.0. The older 
results of the author and his associates, while scattered, agree with the line 
within 12 percent to ug t/Zo=18; the more recent results begin to rise above 
the line at u/ ¢/Z»9=0.5 and are in fair agreement with Hall's data on a 1-inch 
screen. Unfortunately, data at large values of x/M could not be obtained in 
the recent experiments. 

Thus, even when attention is confined to the behavior of the isotropic 
turbulent field, there remain discrepancies in the experimental data such that 
no definite conclusions can be drawn as to the merits of any theory. Further 
experiments are required under carefully controlled conditions in an air 
stream of low turbulence over a wide range of values of x/M and with due 
regard to the various systematic errors that may be present. These experi- 
ments would be of the greatest value if the scale were also measured. 

18. The length \. Relation between \ and L. The general expression for 
the mean rate of dissipation in the flow of a viscous fluid is: 


; f ou\? dv\?2 dw\? /dav ou\? 
romana 
t \ax ay dz dx dy 


(18.1) 


aw  dv\? du dw? 
+(—+—)+(—+— 
oy 02 Oz Ox 


For isotropic turbulence this becomes: 


where up is the viscosity. 
W /du\? ou\? dv Ou 
— = 6(—) +6{—}+6 — — (18.2) 
Ml Ox dy Ox dy 


which, from the relations given in section 14, reduces to: 


W = —7.5p'u?gg = 7.5u(du/dy)*. (18.3) 
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me ia 
go = — 2 lim (18.4) 
r—0 r? 


and has the dimensions of the reciprocal of the square of a length. Let 


But go is defined by: 





go = — 2/r (18.5) 


the factor 2 being introduced to conform to Taylor’s definition of \. Then 
W = 15uu?/d?2. (18 .6) 

The length A may be interpreted in several ways. Equation (18.6) may be 
considered a definition, \ being regarded roughly as a measure of the diame- 
ters of the smallest eddies which are responsible for the dissipation of energy. 
Or, since 1/A?=lim,» (1 —g)/r? =limy—o (1 — R,) Y?, A? is a measure of the ra- 
dius of curvature of the R, curve at Y=0. Or, if a parabola is drawn tangent 
to the R, curve at Y=0, this parabola cuts the axis at the point Y=X. 

Since W= —(3/2)p(du?/dt), the decay law may be written: 


du?/dt = — 10vu?/d?. (18.7) 


This result can also be derived directly from equation (16.5) as shown by 


von Karman, 
By comparing this expression for the decay law with that previously given 


(equation 17.2), namely 


du?/dt = — Au’*/L (18.8) 
it is seen that 
Au’/L = 10v/d? (18.9) 
or, since u’/L/v= No 
?/L? = 10/A No. (18.10) 


Introducing the experimental value of A, 
N/L = 6.97/\/Np. (18.11) 


A similar relation holds for \/L, where L; is the longitudinal scale. If 
the Reynolds Number is formed from L,, the numerical constant is approxi- 
mately 4.93. 

During the decay of self-preserving turbulence No is constant and X is 
proportional to LZ but the constant of proportionality varies inversely as ~/ No 
for turbulent fields of different Reynolds Number. 

Although it cannot be expected on physical grounds that these relations 
hold at very low values of No, there is no experimental evidence of any de- 
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parture from equations (18.7) and (18.8) for values of No as low as 10. There 
seems to be no difficulty in drawing correlation curves for which J is greater 
than L, but no such experimental curves have been measured. However, in 
an example quoted by Taylor,” A/L is as great as 0.86. 

19. The spectrum of turbulence, relation between spectrum and correla- 
tion. The description of turbulence in terms of intensity and scale resembles 
the description of the molecular motion of a gas by temperature and mean 
free path. A more detailed picture can be obtained by considering the dis- 
tribution of energy among eddies of different sizes, or more conveniently the 
distribution of energy with frequency. Just as a beam of white light may be 
separated into a spectrum by the action of a prism or grating, the electric 
current produced by a hot wire anemometer subjected to the speed fluctua- 
tions may be analyzed by means of electric filters into a spectrum. 

The mean value of u? may be regarded as made up of a sum of contribu- 
tions u’F(n)dn, where F(n) is the contribution from frequencies between 
and n+dn and J F(n)dn =1. The curve of F(m) plotted against n is the spec- 
trum curve. According to the proof given by Rayleigh and quoted by Taylor” 


F(n) = 28° lim (I; + 12)/T (19.1) 


Tx 


where T is a long time and 
T 
I, = (1/m) f u cos 2rnt dt 
0 . 
: - (19.2) 
I, = (1/n) f usin 2rnt dl. 
0 


When the fluctuations are superposed on a stream of mean velocity U and 
are very small in comparison with U, the changes in wu at a fixed point may be 
regarded as due to the passage of a fixed turbulent pattern over the point, 
i.e., it may be assumed that 


u = o(t) = o(x/U) (19.3) 


where x is measured upstream at time ¢=0 from the fixed point. The correla- 
tion R, between the fluctuations at the times ¢ and /+x/U is defined by 


7 oot + «/U) 


s (19.4) 
u? 
It can be shown”? that 
f ‘oot + x/U)dt = an’ f (1; + I:) cos (24nx/U)dn (19.5) 
/ «2 0 





20 Taylor, G. I., The spectrum of turbulence, Proc. Roy. Soc. London, Ser. A, 164, 476 
(1938). 
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or, substituting for J]+J} its value in terms of F(n), 


R, = f F(n) cos (2anx/U)dn (19.6) 
0 
and 


F(n) 


(a/v) f R, cos (2anx/U)dx (19.7) 


In other words, the correlation coefficient R, and UF(n)/+/8z are Fourier 
transforms. If either is measured, the other may be computed. R,; is the func- 
tion denoted by f in section 13. The length A, which was defined in terms of the 
function g or R,, is related to R, by the equation: 

1/? = 2 lim (1 — R,)/x?. (19.8) 
z—0 
When # and x are small, cos (27nx/U) in (19.6) may be approximated by 


1 — 21*x?n?/U?. Hence 


1/r\? = (492, us f nk (n)dn. (19.9) 
0 
If the turbulence is self-preserving, the shape of the correlation curve is a 
function of the Reynolds Number of the turbulence. Hence the spectrum 
curve is also a function of the Reynolds Number of the turbulence. Introduc- 
ing the longitudinal scale L, (L.=. R.dx) in equation (19.9), 


i * (nL,\* UF (n) nL, 
soe int f ( ) ——_—— i(~=) (19.10) 
? o \U Lz U 


and in equation (19.7), 


UF (n) 2 2xrnL, Xx x ; 
- -1f R,cos ——— = a( ) (19.11) 
' 0 U a he 


both of which are expressed in terms of the non-dimensional variables 
UF(x)/Lz, nLz/U, x/Lz, }/Lz, and R,. The mean speed U enters only in 
fixing the frequency scale. 

Typical spectrum curves determined experimentally*':” are shown in 
Fig. 10. Studies of the relation between the spectrum and the correlation 
curve have been given by Taylor.”° 

From equation (19.10) it may be inferred that if the curve of UF(n)/L,z vs. 


21 Simmons, L. F. G., and Salter, C., An experimental determination of the spectrum of 


turbulence, Proc. Roy. Soc. London Ser. A, 165, 73 (1938). 
22 Dryden, H. L., Turbulence investigations at the National Bureau of Standards, Proc. Fifth 


Inter. Congr. Appl. Mech., Cambridge, Mass., 1938, p. 362. 
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nL,/U is independent of U, L./X should also be independent of U, which is 
contrary to the known dependence of L,/A on the Reynolds Number of the 
turbulence. Equation (19.9) shows that the value of A is determined largely 
by the values of F(m) at large values of nm. The NPL measurements in Fig. 10 
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Fic. 10. Comparison of National Bureau of Standards and National Physical Laboratory 
measurements of the spectrum of turbulence, plotted non-dimensionally. 
At left, NBS values 40(-) and 160(+) inches behind 1-inch mesh screen at 40 ft/sec. 
At right, NPL values of (F() from Table II of reference 21, L, from reference 20) 82 inches 
behind 3-inch mesh screen at 15(-), 20(X), 25(+), 30(A), and 35({D) ft/sec. 
The reference curve in each case is the curve 
UF(n) + 
i: 4r?n? lL? 
1+ a" 
where U is the mean speed, L, is the integral Sc R.dx, R, is the correlation between the fluctua- 
tions at two points separated by the distance x in the direction of flow, n is the frequency, and 
F(n) is the fraction of the total energy of the turbulence arising from frequencies between m and 





n+dn. 


show clearly this dependence of the spectrum curve on U at high frequencies. 
When the Reynolds Number of the turbulence is large, \/Zz becomes 
small. Experimental measurements show that both R, and Ry, curves ap- 
proach exponential curves. From integration of equation (13.9) it follows that 
2L=L,and equation (19.11) for the corresponding spectrum curve becomes: 
UF (n) 4 
Ls 1 + 49°n?L,/U? 





(19.12) 
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This is the reference curve drawn in Fig. 10. As U decreases, \ increases, and 
the departures at large values of nL,/U becomes greater. The changes in the 
total energy of the fluctuations associated with these changes in the spectrum 
at high frequencies are extremely small. 

Adopting this expression for the spectrum curve, it is possible to compute 
the effect of varying the cut-off frequency of the measuring equipment on the 
measured value of the energy of the fluctuations. If the equipment passes 








03|—\-—_|__— 














04 }——___—+4 Sa ee 4 } : 








02 }-—_/—-~ 





ee ee 











t¢) ol 02 03 04 os 06 o7 os 03 10 


nLy/U 


RATIO OF OBSERVED TO TOTAL ENERGY 


Fic. 11. Effect of cut-off frequencies of apparatus on observed energy of turbulence for 


spectrum given by reference curve of Fig. 9. 
no is the lower cut-off frequency, ms the upper cut-off frequency, L,z the longitudinal scale, 


U the mean speed. 


high frequencies but cuts off sharply at a lower frequency mo, the measured 


total energy is 


—« f* 4(L,/U)dn 4 2anoLs\. — 
set f - ={1- tan~ —)3pu? - (39:43) 
noLe/U 1 + 4a?n?L2/U? 2r U 


The ratio of the observed to the actual total energy is shown in Fig. 11 for 
various values of moL,/U. 

Similarly, if the equipment passes low frequencies but cuts off sharply at a 
higher frequency ,, the measured total energy is (4/27) tan! 27n,L2/U(3pu’). 
The ratio of the observed to the actual total energy for high frequency cut- 
off is also shown in Fig. 10. 

The fact that the correlation and spectrum curves are of the exponential 
type has been interpreted” as meaning that turbulence is a generalized chance 
phenomenon, as nearly chance as a continuous curve can be and retain its 


continuity. 


1943] A REVIEW OF THE STATISTICAL THEORY OF TURBULENCE 37 


20. Fluctuating pressure gradients. In theories of the effect of turbulence 
on transition in boundary layers, it is desired to know the value of the root- 
mean-square pressure gradients, i.e., (@p/dx)*, (0p/dy)*, and (dp/dz)*. Taylor 
has shown® that 





V (ap/dx)? = 2y/2 pu?/d. (20.1) 
Combining this with the relation (18.11), i.e., A/L =6.97/V/(42)!2L/p 





/) ap\?_ 2% oe o/' 5 i (20.2) 
/ (<*) — 6.97L Y i 


The quantities 1/z and L occur in this expression in the combination 
[(\/u2)/L"5|5/2. The ratio (./u?/U)(D/L)"5, where U is the mean speed and 
D the reference dimension of a body under study is known as the Taylor 





turbulence parameter. 

21. The diffusive character of turbulence. An early experimental distinc- 
tion between turbulent and non-turbulent flow was based on the observation 
that a filament of dye introduced into a turbulent fluid stream is rapidly 
diffused over the entire cross section of the stream whereas in a non-turbulent 
flow the filament retains its identity although it may show some waviness. 
It has been pointed out in section 7 that the effect of the turbulent fluctua- 
tions on the mean motion is the introduction of eddy stresses associated with 
the transfer of momentum by the diffusion of fluid particles. Von Karman* 
has given a useful account of the mechanism of the diffusion of discrete par- 
ticles and its effect in producing a shearing stress. A theory of diffusion by 
continuous movements has been developed by Taylor." The process of diffu- 
sion has been found helpful in the experimental study of the statistical prop- 
erties of turbulence. 

22. Diffusion by continuous movements. Consider in a uniform isotropic 
turbulent field the displacement X and velocity u parallel to the arbitrarily 
selected x axis. The intensity +/22 is constant, the field being assumed 
uniform. Let u, and uy be the values of u at times ¢ and ¢’ respectively. 
Consider the definite integral /(u,u dt’. Introducing the correlation coefh- 
cient R,,., between u, and uy, remembering that uv? is constant, 


t t 
fi wacar = Rea, (22.1) 
0 0 


‘'—t=T and place R:,.=Rr. Since Rr is an even function of T, (22.1) ; 


Let ¢’-—t= 
_ s t 
uuydl’ = 7 f RrdT. (22.2) 
J o 0 


may-be written: 
, Th. von, Turbulence, Jour. Roy. Aeron. Soc. 41, 1109 (1937). 
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But 


t t 
f uu, dt) = uf uy dt! = u,X 
0 0 


ux. (22.3) 


Hence 
af RrdT = uX = (1/2)dX*/dt. (22 .4) 
0 
When the time ¢ is so small that Rr approximates unity, equation (22.4) 


becomes: 


(1/2)dX?/dt = ut 
or 


VX? = Vu. (22.5) 


If Rr is equal to zero for all times greater than some time 7 
—_—— To 
uX = af RrdT = constant. (22 .6) 
0 


Define a length J; by the relation: 


To 
L=y ef RrdT (22.7) 
whence ‘ 
ly/ut = uX = (1/2)dX*/dt (22.8) 
and 
X? = hyy/ wt. (22.9) 


If Rr=e-7/T, |, =+/u?T> and the solution of (22.4) yields: 
X? = 2u?T[t — To(1 — e-7/™)]. (22.10) 


Equation (22.10) reduces to (22.5) when ¢ is small compared to 7» and to 
(22.9) when ¢ is large compared to T». 

The diffusion in a uniform field is accordingly completely determined by 
the correlation function Rr. 

23. Diffusion in isotropic turbulence. The foregoing theory is directly ap- 
plicable to diffusion in a uniform isotropic field. However no general state- 
ment can be made as to the relation between the length /; and the scale LZ 
defined in terms of the correlation coefficient R, in section 9. For the turbu- 
lence behind a grid or honeycomb, Taylor found from an analysis of the avail- 
able experimental results that LZ was approximately twice h. 
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The essential features of diffusion in isotropic turbulence expressed in 
equations (22.5) and (22.9) may be summarized as follows: 

1. For time intervals which are small in comparison with the ratio of J; 
to \/u*, the diffusing quantity spreads at a uniform rate proportional to the 
intensity \/, and the rate is not dependent on the length /,. 

2. For time intervals which are large in comparison with the ratio of J; 
to x2, the diffusing quantity N spreads in accordance with the usual diffu- 


sion equation 


ON ON ON ON 0 ON 0 ON fe) ON 

LI a OS ek 

ol Ox oy Oz Ox Ox Oy Oy Oz 02 
with a coefficient of diffusion D equal to hy/u?, where |, is a length defined 
by f, RrdT. 

3. For intermediate time intervals, the diffusion is dependent on the func- 
tion Rr which represents the correlation between the speed of a particle at 
any instant and the speed of the same particle after a time interval 7. 

Consider the diffusion of heat from a hot wire placed in a uniform field of 
isotropic turbulence in a fluid stream of mean speed U. Observations of the 
lateral spread of the thermal wake at a distant x downstream may be used to 
compute the root-mean-square lateral displacement \/Y? of the heated par- 
ticles during a time interval t=x/U. 

It is convenient to characterize the spread by the angle subtended at the 
source by the two positions where the temperature rise is half that at the 
center of the wake. There is a lateral spread of heat produced by the ordinary 
molecular conduction corresponding to an angle ao in degrees of 190.8./ k/pc Ux 
where &, p, and ¢ are thermal conductivity, density, and specific heat (at con- 
stant pressure) of the fluid. It may be shown that the total subtended angle a 
is related to the angle a; produced by turbulent diffusion and qo, as follows: 








2 9 


a =a, + @ (23.1) 


The temperature distribution in the wake follows an “error” curve as does 
the amplitude of the turbulent velocity fluctuations, so that the lateral dis- 
placement Y also has the same Gaussian frequency distribution. The value 
of the lateral spread at which the ordinate is half the maximum is 2.354y Y? 
for this distribution. Hence, expressing @; in degrees, 


a, = 134.7\/¥?2/x (23.2) 
whence from (22.5) for small values of x, 
a, = 134.7/9?/U (23.3) 


where v is written in place of u in (22.5) since the diffusion in the v direction 


is being studied. 
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Thus an experiment on thermal diffusion provides a method of measur- 
ing \/v?. The method was used by Schubauer™ who showed that a, was inde- 
pendent of speed over the range 10 to 50 ft/sec and also independent of x 
over the range 1/2 to 6 inches. 

From measurements at large values of x, it is theoretically possible to 
compute the correlation curve, Rr, vs. JT. In any actual experiment, how- 
ever, the intensity of the turbulence will decrease with x to an extent that 
must be considered. As discussed by Taylor,* Rr may then be considered a 
function of n=Jov vdT =f" (x y?/U)dx. The equation analogous to (22.4) for 
v and Y becomes: 


{ 
\ 


Nz 
(1/2)(U/\/v8) (d¥2/dx) -{ R,dn. (23.4) 
0 


The correlation is given by the expression 


df U 4ayY? 
R, =- ( - -) (23.5) 
dn 2 o as 


and thus involves a double differentiation of experimental curves, a process 
which is usually not very accurate. 

24. Statistical theory of non-isotropic turbulence. In non-isotropic turbu- 
lence the description of the state of the turbulence becomes much more com- 
plex. The eddy shearing stresses do not vanish and the eddy normal stresses 
are not necessarily equal. Six quantities instead of one are required to specify 
the intensity. Similarly the correlation tensor cannot be expressed in terms 
of a single scalar function. In general six scalar functions are required. No 
theoretical investigation using these twelve functions has yet been carried out. 

The exploration of this field is still in its earliest stages. Von Karman?:'® 
has given some discussion of energy transport and dissipation and vorticity 
transport, neglecting the triple correlations, and he has also presented a more 
detailed discussion of two-dimensional flow with constant shearing stress 
(Couette’s problem). The advance of the theory is definitely handicapped by 
the absence of reliable experimental data on the twelve functions required to 
describe the state of turbulence. 

25. Diffusion in non-isotropic turbulence. The only theoretical approach 
at present available for estimating the diffusion in non-isotropic turbulence 
is to consider the process as approximately equivalent to diffusion in isotropic 
turbulence of intensity equal to \/v? and scale 1’, where v is the component 
in the direction in which the diffusion is studied and 1’ is the length defined 
by an equation analogous to (22.7), namely, 


24 Schubauer, G. B., A turbulence indicator utilizing the diffusion of heat, Tech. Rept. 
Nat. Adv. Comm. Aeron., No. 524 (1935). 
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l= vi f Rrdt. (25.1) 
0 


In most experiments the length /’ is not measured. Prandtl defined a mix- 
ing length / in terms of the shearing stress 7 by the relation: 


dU|dU 
dy | dy 


2 


T= 


(25 .2) 








This relation may be interpreted as an equation governing the diffusion of 
momentum with a coefficient of diffusion equal to ?|dU/dy|. Prandtl in fact 
assumed 1/2? proportional to 1|dU/dy| and incorporated the factors of pro- 
portionality in the length /. It is obvious that 

?| dU/dy| = I'v. (25.3) 

The length / can be obtained experimentally if the distributions of velocity 
and shearing stress are known, and, if \/v* is also measured, 1’ may be com- 
puted. Sherwood and Woertz* have made an experimental study of these 
relationships. 

Taylor®® pointed out that fluctuating pressure gradients influence the 
transfer of momentum and suggested that the vorticity be taken as the prop- 
erty undergoing diffusion. The result was the well known vortex transport 
theory. 

Both theories imply diffusion for a time interval long compared to 1'y/2?. 
When diffusion is studied near the source, experiment shows” a behaviour like 
that discussed in section 23. The spread is nearly linear with x, although un- 
symmetical in this case. It is probable that the unsymmetrical character can- 
not be explained on the basis of a single scalar diffusion coefficient. 

26. Correlation in turbulent flow through a pipe. Taylor** has shown that 
the correlation between the component of velocity at a fixed point and that 
at a variable point in the same cross section must be negative for some posi- 
tions of the variable point, if the applied pressure difference between the ends 
of the pipe is constant and the fluid may be considered incompressible. Sup- 
pose the mean velocity is U and the correlation R has been measured between 


25 Sherwood, T. K., and Woertz, B. B., Mass transfer between phases, role of eddy diffusion, 
Ind. Eng. Chem. 31, 1034 (1939). 

26 Taylor, G. I., Transport of vorticity and heat through fluids in turbulent motion, Proc. Roy. 
Soc. London Ser. A, 135, 685 (1932). 

27 Skramstad, H. K., and Schubauer, G. B., The application of thermal diffusion to the 
study of turbulent air flow, Phys. Rev., 53, 927 (1938). Abstract only. Full paper not published. 
A few additional details are given in Dryden, Hugh L., Turbulence and diffusion, Ind. Eng. 
Chem. 31, 416 (1939). 

28 Taylor, G. I., Correlation measurements in a turbulent flow through a pipe, Proc. Roy. 
Soc. London Ser. A, 157, 537 (1936). 
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the component 1 of the fluctuations at a fixed point P and wz at a variable 
point Q in the same cross section. Since the mean flow is constant, 


fw + us)dydz = f Udydz = constant (26.1) 


where the integration is taken over the cross section. At any instant, 


f ued ydz = 0. (26.2) 


Multiplying by m, which is constant for this integration and may be 
placed under the integral sign, 


ff rnsdyas = 0. (26.3) 


Since (26.3) is true for any instant, it is true for the integral over a time 


. 
a/7) f | fsa | dt = 0. (26.4) 
0 


Changing the order of the integration and remembering that (1/T) J uyuedt 


interval 7. Hence 


= Uy Ue 
f Uyu2d ydz = 0. (26.5) 
Introducing the correlation R, 
ff Rus uz dydz = 0. (26 .6) 
But uj is constant with respect to the integration and accordingly 
ff Rusayas = 0. (26.7) 


Since uz is positive, R must be negative for some positions of Q. 
For a circular pipe (26.7) becomes: 


f u; Rrdr = 0 (26.8) 
0 


where u/ is the value of \/u? at radius r and a is the radius of the pipe. 
This relation was experimentally verified in experiments made by Sim- 
mons with the fixed point at the center of the pipe. 
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ON THE MOTION OF A PENDULUM IN A 
TURBULENT FLUID* 


BY 


cc. c tm 
Guggenheim Laboratory, California Institute of Technology 


1. Introduction. In a recent paper, Schumann! has investigated the mo- 
tion of a damped pendulum in a turbulent fluid by considering the effect of 
the fluid as a continuous fluctuating force. He first considers a damped 
pendulum which is bombarded by pellets of equal mass mh at equal intervals 
h of the time, and then treats the case of continuous fluctuations by a limiting 
process. For this latter case, which “must be regarded as being of far more 
practical importance,” Schumann obtains the very interesting result :? 


1 * . 
a | f {R(x-+£)+R(x—£)}e* sin (Bx-++-y)dx 





g 
+f { R(E—x) —R(x—£)} > sin (oxtydx], (1.1) 


The notation is as follows: 
r(£)=correlation function of the displacements of the pendulum at two 


instants separated by a time interval §; 
R()=limiting correlation function of the velocities of the impinging 


pellets; 

\=l+m/M, M being the mass of the pendulum and 7 its damping 
factor; 

8?=0?+-/?—)*, 2x/a being the (damped) period of the pendulum; 
sin y=8/B1; 


B= B+ =08 +P. 

The analysis used by Schumann is very elegant, but somewhat lengthy. 
In this article, we shall study the problem from another point of view, and 
give an alternative derivation of (1.1). This derivation, though unable to 
cover the case of discontinuous impacts, seems to show the nature of that 
relation much more clearly. 

2. Damped pendulum under the action of a fluctuating force. We shall 
now investigate the correlation of displacement of a damped pendulum in 
relation to that of the exciting force. The notation used in this section should 
first be regarded as having different (though analogous) interpretations from 

* Received Dec. 4, 1942. 

1 Schumann, T. E. W., Phil. Mag. (7), 33, 138-150 (1942). 

* Loc. cit., eq. (57), p. 146. 
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those used in §1. The identification of the two systems of notation will be 


made in §3. 
Consider the equation of motion of the pendulum 


d*y dy s , 
—— + 2— + (6° + d*)y = P, (2.1) 
dt? dt 

where y is the displacement, A is the damping factor and 27/8 is the (damped) 
period of the pendulum, P is the exciting force per unit mass of the pendulum, 
and ¢ is the time. If the force P is quasi-periodic, and is given by the real 
part of 

P = >) Ane, (2.2) 


—— 
where w, and A, are real and complex constants respectively, the steady- 
state displacement of the frequency w,/27 is given by the real part of 
a, 
Ye = G,6%%s!, Bg Serres + 
(B* + AXA* — w2) + 2tro,, 
Thus, we have 





= —— (2.4) 
(B? + A? — wi)? + 470i, 


| tat” 
This is the relation between the spectrum of the displacement and that of 
the force in the case of discrete spectra. It is not difficult to generalize this 
result to the case of continuous spectra by the considerations of generalized 
harmonic analysis. We have then 


f(o) = A ———_—_——__ 


where f(w) and F(w) are the spectra of the displacement and the force re- 


f f(w)dw = 1, f Flw)dw = 1 (2.6) 


and A is a constant of normalization, 


spectively, 


These are the well-known relations in the phenomena of resonance. 
The correlations r(€) and R(&) of the displacement and the force respec- 
tively stand in Fourier transform relations to the spectra‘ (apart from con- 


ventional numerical factors): 


3 Wiener, N., “The Fourier Integral” (Cambridge, 1933), p. 150. 
* Wiener, N., loc. cit., eq. (21.21), p. 161, and discussions on p. 163. To be exact, we 
should follow Wiener in calling f(w) and F(w) the spectral densities. 
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2° « 
{(w) = =f r(x) cos wx dx, r(é) -{ f(w) cos wi dw; (2.8) 
Tv 0 


0 


2 -) o 
F(w) = —{ R(x) cos wx dx, R(é) = f F(w) cos wt dw. (2.9) 
T 0 0 


From (2.5), (2.8), and (2.9), we have at once 


2A - cos wkdw - 
r(t) = f - f R(x) cos wx dx. (2.10) 
~ 0 (Bp? + A? — w?)? + 42? 0 








It is not difficult to justify a change of the order of integration, since the cor- 
relation functions are expected to go to zero at infinity sufficiently rapidly. 
The above relation then becomes 


2A " ° COS wk COS wx 
r(t) = —f R(a)dx f dw. (2.11) 
wT Jo o (8B? + dA? — w?)? + 4A%? 





Since we have 
cos wt cos wx = 3} cos w(t + x) + cos w(E — x)}, (2.12) 
we can evaluate the integral with respect to & in (2.11), if we know 


cos wt 


I(t) -f dw 








(2.13) 


This integral is relatively easy to evaluate. We write 


1 2 iwt 
I(t) = =f : de. (2.14) 
2 J_. (8? +? — w)? + 4r%0? 





and consider the corresponding contour integral in the complex w-plane, the 
contour being the usual one composed of the real axis and a semi-circle at 
infinity. The circle is taken in the upper half-plane if ¢>0, and in the lower 
half-plane if t<0. Since (2.13) shows that I(#) is an even function, we shall 
carry out the calculations for ¢ >0 alone. 

There is no difficulty in showing that the integral over the semi-circle 
goes to zero. For, when the imaginary part of wt is positive, | et! is bounded, 
and | (82-+A?—w?)?+ 4\%w?| = O(|w| ) for large values of | | . The evaluation 
of (2.14) then reduces to the calculation of the residues of the integrand at the 
two simple poles +8+iA (A>0) inside the contour. The result can be easily 
verified to be 
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Te 
I(t) = —— sin (Bi , ie 
am Bt + y) (2.15) 


kb 


where 


Pp ; 
(B2 vi x22 ° (2.16) 


With the help of (2.12), (2.13), and (2.15), the equation (2.11) becomes 





sin y = 


A - 
r(é) = =| f ) R(x + &) + R(x — £) fe sin (Bx + y)dx 


+f )R(E — x) — R(x — £)fe* sin (Bx + vax]. (2.17) 


This is Schumann’s relation (1.1), if the constants can be identified. There 
is no difficulty with the normalization coefficient. We have 
46d Tite 
, = 2 R(x)e* sin (Bx + y)dx 
, 0 


e F(w)dw 
= spr f are aaa aera et 
0 (B2 + 2 — w?)? + 4%? 
on putting €=0 in (2.17) and recalling (2.7). The second relation is a by- 
product of our investigation. The limiting case A—0 reduces to the well- 


(2.18) 


known relation (2.9). 

Referring to (2.9), we see that R(£) is an even function of &, so that the 
second integral in (2.17) may be dropped. 

From the derivation, we see that (2.17) is nothing but the Fourier trans- 
form of the well-known resonance relation (2.5). 

3. Identification of the results. We shall now identify Schumann’s result 
with ours by showing that in the limiting case, his pendulum has an effective 
damping factor \ instead of /, and that his correlation of velocity of the pellets 
becomes the correlation of force. The equations of motion as given by Schu- 


mann are® 


M — mh 2mh ; 
; (3.1) 


SS 2 eee 


— Uy + Ur, 
M + mh M + mh 


at the rth impact at the instant ¢=rh, and 
d*y dy 
— + 2/— + (a? + /*)y = 0, 
dt? dt 


(3.2) 


It seems that there are some misprints in the original paper. 


bh 


§ Loc. cit., eqs. (1), | 
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between successive impacts. In (3.1), wf and u, denote respectively the veloci- 
ties of the mass M just before and just after the rth impact, and 2, is the veloc- 
ity of the rth pellet. 

In the limiting case, the discrete impacts become a continuous force given 
by 


Ur — Us 








P,; = lim M = — 2mu + 2m2, (3.3) 


h—0 h 


where u and v are the limiting values for u, and v, at the instant ¢=rh. Thus u 
is evidently the velocity of the pendulum dy/dt, and v is the velocity of the 
infinitesimal impinging pellet in the interval (t, ¢+-d?). The equation of motion 


(3.2) then becomes 


d*y -.. 1 dy 
— + 27I— + (e°® + P)y= _(- 2m —+ 2m2). (3.4) 
dt? dt M dt 


We see that this limiting case carries an inherent damping factor m/M in 
addition to the damping factor /. The term 2mv/M is evidently the exciting 
force per unit mass, since m is the impinging mass per second, the factor 2 
corresponding to the fact that when the impinging mass is infinitesimal com- 
pared with the colliding mass, the former rebounds with the colliding speed,— 
a fact often used in the kinetic theory of gases. 

The identification of Schumann's result with that given in §2 is therefore 
complete. 

4. Discussion. In view of the above derivation, we must be a little careful 
in applying Schumann's relation (1.1) to the study of the motion of a pendu- 
lum in a turbulent fluid. The correlation function R obtained (by suitable 
processes) according to that relation from the correlation function r of the 
displacements is that for the hydrodynamic force (the part corresponding to 
an extra damping being removed). In a turbulent fluid, the connection be- 
tween the velocity fluctuations of the fluid before the introduction of a pendu- 
lum and the fluctuating force acting upon the pendulum after its introduction 
is not at once evident. A preliminary careful investigation seems to be neces- 
sary before the method suggested by Schumann could be used with advantage. 

It is also clear that the spectrum of velocity fluctuations of the pendulum 
will be proportional to w*f(w). It is then clear from (2.8) that the correlation 
function of the velocity fluctuations will be proportional to —r’’(&). 

There is another point which should be mentioned. The spectrum and 
correlation function discussed above refer to those observed at a fixed spatial 
point, if we assume the pendulum to be never far from its position of equi- 
librium. It is not difficult to observe this spectrum with a hot-wire anemome- 
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ter, as has been done by several observers.* The Fourier transform of this 
spectrum will then give the correlation of the fluctuations at successive in- 
stants at the same point in space. It is not clear, however, how (as Schumann 
suggested) the pendulum can be used to observe the correlation function de- 
fined with respect to the same material point,—the quantity introduced by 
Taylor’ in the Lagrangian description of turbulence for the study of the 
phenomenon of turbulent diffusion.® 

In conclusion, the author wishes to express his sincere thanks to Professor 
Theodore von Karman for suggesting the problem to him and for his invalu- 


able suggestions. 





6 For example, L. F. G. Simmons and C. Salter, Proc. Roy. Soc. Ser. A, 165, 73-89 (1938); 
H. L. Dryden, Proceedings of the Fifth International Congress of Applied Mechanics (Cam- 
bridge, Mass., 1938), pp. 362-367; H. Motzfeld, Zeits. f. angew. Math. u. Mech. 18, 362-365 
(1938). 

7 Taylor, G. I., Proc. Lond. Math. Soc., (2), 20, 196-212, (1921). Clearer statements re- 
garding this point are made in his paper of 1935, Proc. Roy. Soc., A, 151, 421-478 (1935). 

8 A discussion of the three types of correlations that may be defined in the study of iso- 
tropic turbulence has been made by G. Dedebant and P. Wehrle, Comptes Rendus 208, 625 
628 (1939). 
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ON PLANE RIGID FRAMES LOADED PERPENDICULARLY 
TO THEIR PLANE* 


BY 


W. PRAGER (Brown University) AnD G. E. HAY (University of Michigan) 


1. Introduction. For purposes of stress analysis, the engineer prefers to 
consider his structures as consisting of plane systems, each of which is subject 
to forces acting in its plane. A typical example is furnished by the conven- 
tional analysis of a parallel chord bridge span in which the side trusses take 
the vertical loads and the top and bottom trusses the transverse loads due to 
wind, etc. In civil engineering this resolution of space systems into plane com- 
ponents is possible in most cases, aifd only very rarely is a structure consid- 
ered as a unit in space. Accordingly, the methods of dealing with space 
structures have not been developed nearly as much as those used in the 
analysis of plane structures. Of course, the general principles of structural 
theory, for instance the principle of virtual work or Castigliano’s principle, 
apply to space structures as well as to plane structures but, as is known from 





Fic. la. Two-sparred wing. 


the case of plane structures, these principles frequently do not offer the most 
convenient approach to the solution of a particular problem. As regards 
special methods, which have been developed so abundantly in the case of 
plane structures, little work has as yet been done in the field of space struc- 
tures. Most of this work is concerned with pin-jointed frameworks. The itera- 
tion procedure of R. V. Southwell’s relaxation method can be applied to space 
structures as well as to plane structures, but efficient direct methods for the 
stress analysis of rigid frames in space are entirely lacking. The present paper, 
intended as a contribution towards the development of such methods, deals 
with the particular case of plane rigid frames carrying loads which act per- 


* Received Dec. 3, 1942. 
1 See R. V. Southwell: Relaxation methods in engineering science, Chap. IV, Oxford 1940, 
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pendicularly to the plane of the frame. An example of this type of structure 
is the monoplane wing of Fig. la, where the spars A—A and B—B are con- 
nected by several main ribs which are fastened rigidly to the spars. Another 
example is the foot ring of an observatory cupola shown in Fig. 1b. 

The method proposed in this 
paper makes extensive use of a 
dual analogy between plane struc- 
tures loaded in their plane and 
plane structures loaded perpen- 
dicularly to their plane. In the 
case of a single straight beam 
this analogy forms the basis of 
the method of conjugate beams? 
which can be considered as a par- 





Fic. 1b. Cupola foot ring. 


ticular case of the present method. 

2. Definitions, notations and sign conventions. This paper is concerned 
with rigid frames consisting of straight or curved members whose axes lie in 
the same plane. This plane is called the structural plane. We will consider only 
frames with members such that every cross section has a principal axis of 
inertia at its centroid lying in the structural plane. Accordingly, when the 
frame is subject to forces acting in the structural plane, the points on the axis 
of any member remain in the structural plane. On the other hand, when the 
frame is loaded perpendicularly to its plane, the displacements of the points 
on the axis of any member are normal to the structural plane. For conciseness, 
the first type of loading will be referred to as plane loading and the second as 
space loading. 

In the case of plane loading the stresses transmitted across any cross section 
of a member of the frame are statically equivalent to the following stress result- 
ants: 1) an axial force which, for the sake of brevity, will be called the pull 
although it may produce either compression or ten- 
sion; 2) a transverse force, called the shear, which , 
acts in the structural plane normal to the axis of 
the member under consideration; 3) a couple, 
called the bending moment, which-also acts in the O x 
structural plane. In the case of space loading the 
stress resultants are: 1) a twisting couple, called 
torque, acting in the plane of the cross section; 
2) a bending couple, called bending moment, whose 
plane is perpendicular to the structural plane as 
well as to the cross-sectional plane; 3) a transverse ee 
force, called shear, whose line of action is normal 
to the structural plane. 








AN | 
W 





Fic, 2. 


2H. M. Westergaard, Deflection of beams by the conjugate beam method, Journal of the 
Western Society of Engineers, 26, 369 (1921). 


1943] RIGID FRAMES LOADED PERPENDICULARLY TO THEIR PLANE 51 


In order to arrive at notations applicable to both types of loading we 
choose on the axis of each member of the frame an origin O and denote by s the 
arc length of the axis measured from this point. A cross section B is then speci- 
fied by giving the corresponding value of s (Fig. 2*). In most cases it will be 
convenient to choose the origin O at one end of the member, in which case s 
will have positive values only. In order to establish appropriate sign conven- 
tions for the loads, displacements, stress resultants and distortions at the 
cross section B, we introduce a rectangular right hand triad with origin at B, 
the x-axis being tangent to the axis of the member at B in the direction of 
increasing s, and the y-axis lying in the structural plane (Fig. 2). 

The loads which the structure carries at B may be forces or couples, either 
concentrated or distributed, or both. The components of the concentrated 
force at B, the distributed force at B, the concentrated couple at B and the 
distributed couple at B we denote by F;, Fy, F:, fz, fy, fa, Cz, Cy, Cay Cay Cy Ce 
respectively, relative to the rectangular triad at B. For example, C, is a con- 
centrated twisting couple, and is positive if its sense is the same as that of 
the 90° rotation necessary to move the y-axis into coincidence with the z-axis. 
Not all the load components thus defined have practical importance; how- 
ever, the analogies which we intend to establish appear more clearly when the 
most general case is considered. 

The force system transmitted across the cross section at B is equivalent 
to a force at the centroid plus a couple. These will be referred to as the stress 
resultants, and we shall denote their components by R,, R,, R., Mz, My, M, 
respectively, relative to the rectangular triad at B. R, is the pull, R, and R, 
the shears parallel and perpendicular to the structural plane, M, the torque 
and M,, M, the bending moments. 

The stress resultants and the loads are connected by the equations of 
equilibrium. If no concentrated forces are applied at the cross section B, the 
equations of equilibrium for a straight structural member are: 


R: + PF = Se M; + Cz _ 0, 
R, + fy = 9, Mj +c,— R,=0, (1’) 
Ri+f.=0, Mitea+R,=0, 


I 


where the dashes denote differentiation with respect to the arc length s, If 
concentrated loads F,, Fy, F,and Cz, Cy, C, are applied at B, we have 


Sit) -Be~ 4h oh B40 ~ BS - 84606 
0, M,(s+e6) — M,(s—e) +C, = 9, (1’’) 
0, Ms+e6—-MAs—6)4+C, 0, 


Ris+6) —R(s-—O+F, 
RAs+e) —RA(s-o+F, 


where € denotes an arbitrarily small length. 


* In Fig. 2 the origin of the system x, y should be marked B. 
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wn 
iS) 


The displacement of the cross section B is specified by the components 
Uz, Uy, Uz Of the translation of the centroid and the components 62, 6,, 0. of 
the rotation of the cross section. 

Finally, the six distortion components, gz, Zy, ge and hz, hy, h,, of a straight 
structural member are defined as follows: 


fs = Us, h, = 6, 
Sy = ty — 4, hy = 4,, (2’) 


ur + 4,, h, = 0, 


II 


Sz 
where the dashes again denote differentiation with respect to the arc length s. 
gz will be called the stretch, g, and g, the slips, hz the twist and h,, h, the bends. 
Two structural members may be connected by a link which permits some 
relative displacement of the end sections of the two members, for instance 
by a hinge permitting a free bend. Such relative displacements can be handled 
as concentrated distortions: 
6.(s + «€) — 0,(s — €), 


G, = u(s + €) — u,(s — ), BH. 
G, = u,(s + €) — u,(s — ©), H, = 0,(s + €) — 0,(s — e), iZ"3 
= u(s+e) — u(s — 6), H, = 0,(s + €) — 0,(s — €), 


II 


G 
where € is again an arbitrarily small length. 

For elastic structural members the stress resultants can be represented as 
the products of the corresponding distortions and stiffness factors: 


\ 


R, = Azrkz, R, = Ayfy, R, = A282, 


z 


M. = 6,4, M, = B,h,, MM, = Bk; 
where a, = EA, a, =GA/ky, a, =GA/k., Bz is the torsional rigidity of the mem- 
ber, 6, = EI,, 8. =EI,, E being Young’s modulus, A the area of the cross sec- 
tion, G the modulus of rigidity, k, and k, constants depending on the shape 
of the cross section, J, and J, the moments of inertia of the cross section with 


(3) 


respect to the axes of y and z. 

3. Analogy between statics of plane loaded and kinematics of space 
loaded frames. For a plane loaded frame the loads F,, f., Cz, cz, Cy, cy and 
the stress resultants R,, Mz, My vanish. Similarly, for a space loaded frame 
the displacements uz, uy, 6, and the distortions Gz, gz, Gy, gy, Hz, hz are zero. 
The remaining equations (1) for the members of the plane loaded frame then 
are seen to correspond exactly to the remaining equations (2) for the members 
of the space loaded frame according to the Table I. 











TABLE | 
—————— — ———————————————— oF ——————_—————_—————————————————— — — ——————————— 
Statics of Plane Loaded Frame Kinematics of Space Loaded Frame 
Pris Pye ee [-H., —H,, —G 
Loads ( , Distortions . : 
fer fur Ge (hs, —hy, —£s 


Stress resultants R,, R,, M, Displacements  9,, Oy, Uz 
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At a point B where two straight members are rigidly fastened to one an- 
other under the angle ¢ (Fig. 3), we have the following relations between the 
stress resultants of the plane loaded frame at the two sides of B: 


R, = R, cos ¢ + R, sin ¢, R, = — R,z sin ¢ + R, cos ¢, M, = M,. 


These relations correspond exactly to the following relations for the displace- 
ment components of the space loaded frame on the two sides of B: 


6, = 6, cos ¢ + 6, sin ¢, 6, = — 0, sin ¢ + 4, cos ¢, il, = tz. (5) 


A curved member may be considered as the limiting case of 
a polygonal arrangement of straight members. Accordingly, 
the correspondence between the stress resultants of the plane x > 
loaded frame and the displacements of the space loaded frame 
remains valid in the case of frames containing curved mem- 
bers. Of course, equations (1) and (2) in their present forms 
do not apply to frames with curved members. Equations ap- - Ys 
plicable to such frames will be developed in a later paper. 

If more than two members join at the same point, the Y 
analogy breaks down. This can be seen from the example in x* 
Fig. 4, where we have for the stress resultants of the plane 


loaded frame Fic. 4. 


R. = R, — R} R, = R, + Re, M, = M, — M3, 


uv? 


and for the displacements of the space loaded from 


6. = 0, = 0* 6, = 0, = — OF, i-=u,= us. 


: y? 
The analogy therefore applies directly only to frames which consist of a 
simple chain of members without branch points. In spite of this fact the 
analogy is very useful even in the case of more 
._0 b complex frames since these may be considered 
as consisting of simple frames to which the 

¢C analogy can be applied. 
In order that the analogy indicated above 
B hold everywhere in the structures, it is neces- 
sary that there can be a certain correspondence 
between the various supports and links in the 























“” 
—_ . . 
3 Ss frame with plane loading and the frame with 
t 3 space loading. For example, the analogy is 
ait a eR iateiaatics J maintained if to a pin support in the frame 


with plane loading there corresponds a simple 
support in the plane with space loading. A 
consideration of certain types of supports and links leads to the results pre- 
sented in Table II. This table is by no means complete. 


Fic. 5. 
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We shall now use this analogy to determine the influence lines for the 
stress resultants of a statically determinate space loaded frame. We consider 
the rectangular frame with three pin supports 1, 2, 3 shown in Fig. 5. Ac- 
cording to a well known principle of the theory of structures, the ordinates 
of the influence line for the bending moment produced at the section B by 
transverse loads F, are precisely the displacements u, produced by a unit 
bend, H,=1, at B. According to the analogy explained above these displace- 
ments can be obtained as the bending moments, M,, produced in the corre- 
sponding plane loaded frame by the load, Fy=—1, at B. This latter frame is 
called the conjugate frame and is the three hinged portal shown in Fig. 6a. 
The bending moments due to the unit load are easily computed and are shown 
in Fig. 6b. The diagram showing the bending moments of the conjugate 
frame is at the same time the influence line for the bending moment produced 
at the section B of the space loaded frame by transverse loads F,,. 





a/3 a/3 














2/3 1/3 


a/3b a/30 
Fic. 6a. Fic. 6b. 


Similarly, the influence line for the torque M, produced at the section B 
of the space loaded frame by transverse loads F, can be obtained as the bend- 
ing moment diagram of the conjugate frame due to the load F,= —1 at B. 
Finally, the influence line for the shear R, produced at the section B of the 
space loaded frame by transverse loads F, is found as the bending moment 
diagram of the conjugate frame due to the couple C,= —1 at B. 

4. Analogy between kinematics of plane loaded and statics of space 
loaded frames. For a plane loaded frame the displacements u,, 62, 0, vanish 
and the distortions G,, gz, Hz, hz, Hy, hy are zero. For a space loaded frame 
the loads F., fz, Fy, fy, Cz, ¢, and the stress resultants Rz, Ry, M, vanish. The 
remaining equations (2) for the members of the plane loaded frame then are 
found to correspond exactly to the remaining equations (1) for the members 
of the space loaded frame according to Table III. 

The relations between the displacements of the plane loaded frame at both 
sides of the angle joint of Fig. 3 are easily seen to correspond to the relations 
between the stress resultants of the space loaded frame. Furthermore, the 
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correspondence between the various types of supports and links shown in 
Table II is valid also in the present case where we are concerned with the 
kinematics of the plane loaded frame and the statics of the space loaded 





frame. 
TABLE III 
Kinematics of Plane Loaded Frame Statics of Space Loaded Frame 
Displacements tz, Uy, 9% Stress resultants M,, M,, R, 
|} -—G,, —G,, -Z. ic. .. F 
i " : Loads Nak 43 ¥ 


Distortions 
—£2, —Zy, —hs 


yy 


(Cz, Cy, fe 
The analogy established in this section can 
be used to find the stress resultants of a stat- 
ically determinate space loaded frame by 
determining the displacements of the con- 
jugate plane loaded frame. Let us consider 
for instance the frame in Fig. 5 carrying a 
? transverse load F,=FP at B. The stress result- 
ants produced by this load can be found as the 
displacements of the conjugate frame produced 
by the bend H,= —P at B. The general trend 
of these displacements is shown in Fig. 7. Rel- 
Fic. 7 ative to the right hand support, we obtain for 
the transverse displacement at B:u,=4a¢g/ 3. 


ap be, ae? 


| 









a/3 
et 


arse 





2a/3 


Relative to the left hand support, we have for the transverse displacement 
at the same point: uy=2a(P—g)/3. By equating these expressions and 
solving for yg, we find that g=P/3. The displacements u,, uy, 0, of the con- 
jugate frame are thus known. They correspond to the stress resultants 


M,, M,, R. of the space loaded frame, the distribution of which is shown in 


Figs. 8 a—c. 



































| | 1 | TTT] Pb/3 -Pa/3 -P/3 
Be if 1 | | L| Lt 
Pass Pa/3 PT A 4 
rP/3 P/3--4 
+ 2P/3 = 
M, M, R, was 


Fic. 8a-c. 


5. Elastic deformations of space loaded frames. Indeterminate space 
loaded frames. When the stress resultants of a space loaded frame have been 
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determined by the method given in the preceding section, the analogy of §3 
can be used in order to find the elastic deformations of the space loaded frame. 
For example, let us consider again the frame in Fig. 5, carrying a transverse 
load F, =P at B. The stress resultants have been determined in the preceding 
section and are shown in Fig. 8. Equations (3) furnish the distortions pro- 
duced by these stress resultants: h, = M,/Bz, hy = M,/By, g. =R./az. According 
to the analogy of §3, the displacements 6@,, 6,, uz, corresponding to these 
distortions can be found as the stress resultants R,, Ry, M, of the conjugate 
frame carrying the loads f,= —hz, fy= —hy, c: = —g:. In most cases the influ- 
ences of the shear R, on the deformations can be neglected, and consequently 
in the conjugate frame the distributed couples c, need not be considered. The 
loads on the conjugate frame then consist of the axial loads shown in Fig. 9a 
and the transverse loads shown in Fig. 9b, it having been assumed that 6, 
and £, have the same values for all members. The reactions can be computed 
easily and are indicated in Figs. 9a and 9b. If now, for instance, we wish to 
determine the deflection u, of the right hand corner C of the space loaded 
frame, we have only to compute the bending moment J/, at the corresponding 
corner of the conjugate frame. We find that 
Pa 

—— [23 — 18n* + 54yn], (6) 
1628, 


“u,=— 


where n =b/a and y =8,/Bz. 
This method of computing elastic deformations enables us to carry out 


the stress analysis of indeterminate space loaded frames. Let us suppose, for 
instance, that the frame in Fig. 5 is given a further simple support at B and 
carries a transverse load F, =Q at C. By establishing the condition for vanish- 
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A Abd qb 
qa qa Vi % 
Fic, 9a. gq=Pa/38z. Fic. 9b.:n=b/a, r= Pa/3£,, 
ra Po ra G 
U, = — (40 + 18n'), U2, = — (32 — 18m"). 
54 54 
: ra a : ra 75 
V, = — (23 + 270°) Vo = — (23 — 27n'). 
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ing deflection at B and applying Maxwell's law of reciprocal deflections in 
the usual way, we find the reaction at B in the form R= —Quc/ug, where 
up and uc denote the deflections which a unit transverse load F,=1 at B 
produces at the points B and C respectively; uc can be obtained from (6) by 
setting P=1; uw, can be computed as the bending moment WM, at the point 
B of the conjugate frame, loaded according to Figs. 9a and 9b. With P=1 
we find that 
a® 
un = —— [32 + 18% + 54yn | 
2438, 


and therefore 
3(23 — 18n* + 54yn) 





2(32 + 18° + 54yn) 


6. The inverse column analogy. The following method of determining the 
stress resultants of indeterminate space loaded frames is patterned after 
the column analogy method of H. Cross.’ In 

Le b order to avoid lengthy computations which 











Le ole. 
, might obscure the essential feature of this new 
; a B , é 
method, we shall consider the simple problem 
of the frame in Fig. 10 carrying a transverse 
a load F,=F at B. 
| We suppose again that the stiffness factors 6, 
and 8, have the same values for all members. 
| The frame then is symmetrical and is loaded 
aot oh. symmetrically. 
Af. ° » 
71M Za According to the analogy of §4, the torque 





M, and the bending moment M, of the frame 
can be obtained as the displacements u, and u, 
of the conjugate frame produced by the bend H,= —P at B. Now the conju- 
gate frame is entirely free. Since the system is symmetrical, it will possess a 
kinematically indeterminate displacement v in the direction of the columns as 
shown in Fig. 11a. This indeterminate displacement in the plane loaded frame 
corresponds to the statically indeterminate torque in the columns of the space 


Fic. 10. 


loaded frame. 
From equations (2) we see that the distortions of the space loaded frame 


are given by h.z=M,/8:, hy=M,/B,, or, by the analogy of §4, by 4z=u./:, 
hy =u,/By, where uz and uy, are the displacements in the conjugate plane 
loaded frame. The analogy of §3 then indicates that the longitudinal and 
transverse loads on the conjugate frame are given by f, = —uz/Bz, fy = —Uy/By. 


No. 215, 1930. 
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The horizontal load is in equilibrium because of symmetry. The vertical load 
is shown in Fig. 11b, and is in equilibrium if 


Pb n 


ee ee 

4 nt+y 

where n =b/a, y =8,/8z. Thus the torque M, in the left hand column is given 
by M,=u,=v. Also, the bending moment M, at B is given by 

Pb 2y+n 


M, = — (v — Pb/2) = — 
v \ /2) 4 vy+n 


o 
v= 


The procedure outlined above is equivalent to the following procedure. 
We suppose that the conjugate frame is embedded in an elastic jelly which 











Rarer Pb/2 eee. See 
wR, | | 

_ ] ] 

-~ { 

WT { 

wi 

: y wf 
Pavel Pale | { 

Kd : 
Fic. 11a. Fic. 11b. n=b/a, r=Pa/2£,. 


offers resistance to the displacement of the frame in such a way that an ele- 
ment ds with displacements u, and uy will meet with a resistance consisting 
of longitudinal and transverse forces of magnitudes u,ds/8, and u,dy/B,, re- 
spectively. We then determine the displacements of this elastically supported 
frame, and the longitudinal displacement uz then gives the torque M, in 
the space loaded frame, while the transverse displacement u, gives the bend- 
ing moment M,. The above method of procedure is clearly seen to be the 
counterpart of the column analogy of H. Cross; it will be called the inverse 
column analogy. 

The inverse column analogy furnishes a simple method of determining 
the influence lines of an indeterminate space loaded frame. For example, let 
us consider the influence line for the bending moment M, which transverse 
loads F, produce at the section B of the symmetric frame in Fig. 10. Following 
the line of approach of §3, we consider the conjugate frame loaded by a 
transverse force F, = —1at B and supported elastically in accordance with the 
inverse column analogy. Since this frame is symmetrical and loaded sym- 
metrically, the displacements are of the type indicated in Fig. 12a. The forces 
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which the elastic support exerts on the frame are shown in Fig. 12b. From 
the condition of equilibrium for the vertical forces, we find that 
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where again n=b/a and y =8,/8:z. The bending moment diagram of this elas- 
tically supported frame (Fig. 12c) is at the same time the desired influence 


line of the space loaded frame. 
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ON THE VIBRATIONS OF A CLAMPED 
PLATE UNDER TENSION* 


BY 
ALEXANDER WEINSTEIN ano WEI ZANG CHIEN 
Department of Applied Mathematics, University of Toronto 


The object of the present paper is the computation of the fundamental 
frequency of a vibrating clamped square plate under uniform tension. It will 
be seen that the method used here reduces our problem for a plate of any 
shape to the membrane problem for the corresponding domain. For this rea- 
son similar numerical results could be obtained for a number of other shapes. 
A similar question has been discussed for a circular plate by W. G. Bickley! 
in connection with the problem of reception of acoustic signals in a condenser 
microphone. The circular plate is an elementary problem from the theoretical 
viewpoint. However, the actual calculations involving Bessel’s functions are 
rather heavy. Bickley was able to give the frequencies only for a small range 
of the tension. 

The frequencies of a square plate cannot be obtained explicitly in terms 
of elementary functions. However, the Rayleigh-Ritz method yields an upper 
bound for these quantities. The result cannot be considered as satisfactory 
since this method does not give us an estimation of the error. Fortunately, 
an increasing sequence of lower bounds can be obtained for all frequencies 
by the application of a variational method already introduced by one of the 
authors in several vibration and buckling problems. Combining these lower 
bounds with the upper bounds obtained by Rayleigh-Ritz, we obtain a narrow 
interval in which our frequencies are located. 

Moreover, it is obvious that for questions like that of microphone recep- 
tion, the lower bounds are the more important data. 

The theory of the new variational method has been developed in several 
papers.? The modifications in the present case are slight. For this reason we 
will omit all theoretical details. The reader can easily reconstruct the proofs 
of the rules which we are following here. 

Let S be the domain of a plate of arbitrary shape, and let C denote its 
boundary. In the numerical applications we shall assume that S is the square 
—n/2<x, ySr/2. 

We denote by: 

2h, the thickness of the plate 

7, the tension 

* Received Dec. 11, 1942. 

1 W. G. Bickley, Phil. Mag. (7), 15, 776-797 (1933). 

2 A. Weinstein, Mémorial des Sciences Mathématiques, No. 88, 1937; A. Weinstein, 
Portugaliae Mathematica 2, 36 (1941). 
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E, Young’s modulus 
g, Poisson’s ratio 
p, density of the material 
D =2Eh*/3(1—o?), the flexural rigidity 
w, the eigenfrequency (number of pulsations in 27 seconds) 
w, the transversal displacement. 
We put r=77/D. Our problem admits an infinite sequence of eigenfrequencies 


W =}, We, - * * , in place of which we shall use the eigenvalues A=, Ax, + + - , 
where 
2hpw? 
D 


The displacements corresponding to these eigenvalues will be denoted by 
- + +, These transverse displacements w satisfy in S the differential 


—_— 


9 
am 


equation: 


W, We, 


AAw — rAw — Aw = 0 (1) 


with the boundary conditions 


w= 0 (2) 
dw/dn = 0 (3) 
on C, 
The equation (1) may be written as follows: 
(A+ a)(A — B)w = 0, (a > 0,8 > 0) (4) 
with 
B-—-a=r, aB = X. (5) 
or 
T L ee T 7 : 
i es Sree ; i ne nae : ¢ 
ax 7 1 + A, B , + rl + (6) 
We sce that we have the identity: 
w=utéaéa in S+C, (7) 
where uw and 4@ are solutions of 
Au + au = 0, (8) 
Aa — pa = 0. (9) 


We have therefore also the following identity: 
Aw = A(u+ a) = Bi —au in S+C. (10) 


The identities (7) and (10) will be useful in the following. 
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It is well known that the eigenvalues of our plate can be defined by minima 
problems, the same as could be used in the Rayleigh-Ritz method. For in- 
stance, the fundamental eigenvalue \ =\,, in which we are interested in this 
paper, is given by the variational problem: 


ow \? Ow \? 
U(w) = ff cawranay a “ff (5) + ( -) dxdy 
Ss Ss Ox oy 


= min = dX, 


H(w) ={f wdxdy = 1 (12) 
Ss 


and with the boundary conditions (2) and (3). 

Let us note that U is not the potential energy of the plate. Nevertheless 
our variational problem gives us the correct differential equation and bound- 
ary conditions. This variational problem will be denoted by P. The higher 
eigenvalues As, As, > + can also be defined by similar variational problems. 
However, we shall not use them in this paper. 

The Euler equation of P is the equation (1). This equation together with 
the boundary conditions (2) and (3) defines a differential eigenvalue problem 
P which admits the solutions w;, we, - - - corresponding to the eigenvalues 
a. re 

In order to obtain an increasing sequence of lower bounds for \, we begin 
by cancelling in the variational problem P the boundary condition dw/dn =0. 
In this way we obtain a new variational problem Po: 





(11) 


with the condition 


: ’ (0) 
U(w) = min =), ; H(w) = 1 (13) 
with the boundary condition w=0. 
The conditions in P» being less restrictive than in P, we have df <\x. 
The Euler equation in Po is the same as in P, namely the equation (1). How- 
ever, the boundary conditions for this equation are 


w= 0 and Aw = 0 on C, (14) 


the last condition being a so-called natural boundary condition, i.e., a condition 
which is automatically satisfied by the minimizing function in Po. The corre- 
sponding differential eigenvalue problem Py is given by (1) and (14). Po ad- 
mits a sequence of eigenvalues A, A, ALS’, - - - the smallest of which is 
identical with the minimum X in Po. In the case of a square plate, Po is 
identical with the problem of the vibrations of a supported plate under ten- 
sion. The problem Pp can be solved, for a domain of any shape, in terms of 
the membrane problem for the same domain, a fact which has been implicitly 
used in the elementary theory of a square supported plate. 
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In order to show this we use the identities (7) and (10) and we obtain 
from (14) at once 


“a= 0, @=0 on C. 


In view of (8) and (9), it follows immediately that #=0 in S+C and that (1) 
and (14) are satisfied by w=u, where u is an eigenfunction of the membrane 
problem Au-+-au =0 in S, u=0 on C. From the eigenvalues a of this problem 
we can compute the eigenvalues A in Py by using the equations (5). 

In order to obtain an increasing sequence of lower bounds for the eigen- 
value \; in P we link Po with P by a chain of intermediate variational prob- 
lems P;, Po, - ++, the solutions of which can be expressed in terms of the 
solutions of Po. In this way our problem for the clamped plate can be reduced 
to the solution of Py which is, as we have seen, equivalent to the problem of 
a vibrating membrane. 

In order to show how this can be done, let 


Pils), pols), -- +, Pm—as(s), Pm(S),° °° (15) 


be an arbitrarily given sequence of functions defined on the boundary C of 
the plate, s being the arc, and f:(s) being positive. The problem P,, 
(m=1, 2, -- - ) is then defined as follows: 


Problem P,,: Find the minimum \®” of U(w) with the condition H(w) =1 
and with the boundary conditions 


w=0onC, (16) 


dw a’ 
pr — ds = 0, 2? ae (17) 
c dn 
The conditions in P,, are more restrictive than those in P,,_; but they are 
less restrictive than the boundary conditions in P. 
We have therefore \ <\M <r?) < - - - Sry. The minimizing function w 
in P satisfies the same Euler equation (1) as in Po (or in P), but the boundary 
conditions are given now by the equations (16), (17) and by 


Aw == a;p1 + ga A + QmPm ON G. (18) 


the last condition being again a natural boundary condition. The constant 
coefficients a;, - - - , @m are unknowns. In order to solve P», we have to de- 
termine the lowest eigenvalue of the differential problem P,, defined by (1), 
(16), (17) and (18). We use again the identities (7) and (10) in a way de- 
scribed in our previous papers. In order to avoid repetitions which would 
considerably increase the length of this paper, we will only formulate the 
rules for the computation. It can be shown that \%” can be computed by the 
following procedure. Denote by u; and #; ({=1,---, m) the solutions of 


the equations 
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Au; + au; = 0 (19) 

Aa; — Ba; = 0 (20) 
with the boundary conditions 

u; = — p(s) (21) 

ui; = pi(s) (22) 


where a@ and £ are considered as parameters. These equations can be solved 
in terms of the solutions of Po. Put w;=u;+a; and compute the quantities 


dw; 
a;j(A) = fo ds; i,j=1,--+,m, (23) 
- 


dn 


where the parameter A is defined in terms of a and B by (5) and (6). Then AY” 
is the smallest root of the determinant equation 


|aii(A)|| = 0; 1j=1,--+,m (24) 
provided that the smallest root is smaller than the second eigenvalue d{ of 
the differential eigenvalue problem Po, defined by the equations (1) and (14). 

The calculation of (23) can be further simplified by introducing a sequence 
of harmonic functions 


Pil x, y), po(x, y), dine Pm-1(%, y), Pm(x, y), aes (25) 
whose boundary values are given as in (15). Then by Green's theorem, (23) 
can be written as follows 


a;(\) = f pi(x, y)(Ba; — au;)dxdy. (26) 


Calculation for a square plate (—1/2Sx, ySm/2): In this case, we take 
cosh (82;-19/2) cosh (a2;~:7/2) 
cosh (27 — 1)2/2 


-!cosh (2i — 1)x cos (2i — 1)y + cosh (2i — 1) y cos (2i — 1)x} (27) 





p(x, y) = 


where 
a1 = V/(2i — 1)? -a, Boa = V/(2i — 1)? +8. (28) 

On the boundary, we have 
{p(t m/2,y) = cos (2i+1)y cosh Be:-14/2 cosh ao;-197/2 


(29) 
Voix, + 2/2) = cos (2i — 1)x cosh Bo;_:97/2 cosh ao;-17/2. 


Then the solutions of the problems (19), (21) and (20), (22) are 
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—cosh 62;-17/2[cos (2i—1) x cosh aa:_1y+ cos (2i— 1) y cosh ay;_ x | 


(30) 
cosh ao;_197/2 [cos (2i—1)x cosh B2;-1y-+cos (2i—1) y cosh Boi—14]. 


= = 
= 8 
°. °. 


Substituting p:(x, y) from (27), u; and a; from (30) into (26), we obtain after 
a little calculation, 


a;; = 4 cosh ao;197/2 cosh a2;-19/2 cosh B2;-19/2 cosh B2;-19/2 (Ai; + Bij) (31) 


where 
2(27-—1)(2i—1)(—1)**? B a 
A;j=A p= : ae ee | l= on +-- | (32) 
[(2é—1)?-+(2j—1)?] L632, ,+(2i—-1)? of, + (2¢—1)*L 
pa : [Boi 17 tanh Boi-iT, 2— a9;- 17 tanh Qoj-17T — 1/2] 33 
(33) 
ie for 1¥*7. 


It should be noted that the roots of |/a;;(A)|| are equal to the roots of ||&;(A)]|, 
where 

y= Ay + Bu. (34) 

The results of our numerical computations are given in Table I below. 


The first and second columns give \ and for the supported plate. The 
next three columns give the smallest root of the determinantal equation (24) 








TABLE [| 
Supported Plate Clamped Plate 
ist eigen- | 2nd eigen- | a (3) (3) Rayleigh- 
value \\” | value ogg Me * » | Ritz method 
5 | 14 50 24.982 | 25.222 | 25.236 | 25.509 
10 24 75 36.639 | 36.845 36.862 | 37.443 
2 34 100 48.084 | 48.253 48.284 | 49.261 
20 | 44 125 59.289 | 59.452 59.491 | 61.008 
30 64 175 81.651 | 81.760 81.809 | 84.372 
50 104 275 | 125.43 | 125.56 125.59 | 130.85 
100 204 525 225.56 225.63 225.65 | 246.58 
1025 443.15 443.24 443.25 477.58 


200 404 








for m=1, 2, 3. Since these roots are smaller than the corresponding second 
eigenvalues dS, they are, according to the general theory, identical with the 
eigenvalues \{”, \?, A! and give therefore an increasing sequence of lower 
bounds for the fundamental eigenvalue \; in P. The corresponding upper 
bounds, obtained by the Rayleigh-Ritz method are tabulated in the last 


column. They have been obtained from the variational problem P by putting 


w = A cos x cos* y + B cos® x cos? y. 
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A comparison with \{® shows us that the error in the values of ), is, for small 
tensions, less than 1.2 per cent and, for great tensions, less than 7 per cent. 
The fact that \{° hardly exceeds \? makes it probable that the lower bounds 
are much closer to the true value of \; than the upper bound given by the 
Rayleigh-Ritz method. 

In figure 1 are plotted curves of fundamental eigenvalues of clamped circu- 
lar plate and square plate against the tension r. The curve I is the values of 
\® for the clamped square plate. The curves II and III are respectively the 
fundamental eigenvalues for a circular plate of equal area and equal circum- 
ference as the given square plate. Both of the latter curves are calculated 


from the Bickley result. 
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Fic. 1. Curve 1: Clamped Square Plate (1/2 2x, y= —2/2) 
Curve 2: Clamped Circular Plate (r = x) 
Curve 3: Clamped Circular Plate (r= 2) 


Remark. Using our lower bound \® for a single value ro of t we can easily 
compute lower bounds for \; for every value of r. This result can be obtained 
by combining our method with an idea of R. V. Southwell.’ In fact, the lowest 
eigenvalue \; =A,;(7T) is given by the minimum of U(w)/H(w) under the con- 


3H. Lamb and R. V. Southwell, Proc. Royal Soc. Ser. A, 99, 272 (1921). 
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ditions (2) and (3). Denoting in (11) the first and second integrals by J(w) 
and D(w) respectively we have 


U/H =J/H +7rD/H =J/H + 1oD/H + (7 — 10)D/H. (35) 


Since D(w)/H(w) is obviously greater than 2 (i.e., greater than the lowest 
eigenvalue of the vibrating membrane) we have for all values rT >To 


Ai(7) > Ai(79) + 2(r — To) ba Ar (re) + 2(r = To). 


Putting r>=5 we give in Table II the values of A (5)+2(r—5). 





TABLE II 

r di(5) +2(7 —5) 

5 25.236 

10 35.236 

15 45 .236 
20 55.236 
30 75.236 
50 115.23 
100 215.23 
200 415.23 


It will be seen that these lower bounds for \; are smaller than the lower 
bounds computed by our method alone. 
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A DIRECT IMAGE ERROR THEORY* 


M. HERZBERGER 
Communication No. 894 from the Kodak Research Laboratories 


1. Ina previous paper! the author proposed a direct approach to the prob- 
lem of geometrical optics. In this paper we shall give a new image error theory, 
to the fifth order, which seems to be more adapted to the practical problems 
than former theories. We are given a rotationally symmetric system. Let us 
choose two Cartesian systems, one in object space and the other in image 
space, such that the x, x’ and y, y’ axes have the same directions and the g, 2’ 
axes coincide with the optical axis of the system. 

A ray is given in object (image) space by the coordinates x, y, (x’, y’) of 
its intersection point with the plane z=0, (s’=0). The optical direction 
cosines (the direction cosines multiplied by the refractive indices m and n’, 
respectively) may be designated by the Greek letters 


Ent = Ve — (+97); 8,9, 0 = vn? — (8% + 9”). 
The fundamental problem of practical optics is to find x’, y’, &’, n’, when 
x, y, & » are given. Because of the rotational symmetry, four functions, 
A, B, C, D, exist such that, 
x’ = Ax + Bi, t = Cx + Dé, 
y’ = Ay + Bn, ni = Cy + Dy. 
where A, B, C, D depend only on the three symmetric functions 1, U2, us of 


our coordinates: 


ui = 3(?+ 9°), ue= att yn, us = 3(6? + 7°). (2) 








(1) 


We found in the previous paper! that, according to the laws of geometrical 
optics, A, B, C, D cannot be arbitrary functions, but must fulfill one finite 
and three differential equations, viz., 


AD — BC =1 


and 


dC oD OA OB 0A dC dA OC 
a(- > 2 c(- ==) + 20 (— —=— —— 
Ou» Ou} Ou> Ou, Ou, Ou» Ou. OU; 
(< oD 0A aD OB aC OB _) 
uo - 


+ 


Ou; OUe Ou, Ou, Ou; OU2 Ou2 OU, 


(— oD 0B aD 
2us ) 


+ 


Ou, Ot Ou. OU, 


* Received Nov. 26, 1943. 
'M. Herzberger, Direct methods in geometrical optics, Trans. Amer. Math. Soc., 53, 218 
(1943) 
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oC aD oC oD OA OB 
4{—— — —)])+ Bit — —- —]} -—- C{—_ - — 
Ous Ous Ouse Ou OUu3 Ous 


OA OB 0A dC 0A dC 
-— DO — —- —]) + 2a —- —- - 
Ous Ou, Ou, OuU3 Ouz OU, 


0A aD 0A aD OB a 0B aC 
+ m(— — - 2 een German eee ae 


Ou, Ous Ou3 OU, Ou, OUs3 Ou3 OU, 


0B aD 0B aD 
ss 20 (—— ee ee | oer —) = 0), 


Ou, OUs Ou3 OU, 


dC oD OA OB 0A OC 0A OC 
B(— — —} — Di — —- —) + 2u;| — — -—- — — 
Ous Ouse Ou3 Ouse Ous Os OuU3 OU2 
0A OD 0A aD 0B dC OB dC 
+ Us 


Ou2 OU3 Ou3z Ouse OuUs OU3 Ou3 Oe 


It is the purpose of this paper to develop from formulae (A) the theory of 
image errors. Developing A, B, C, D into a series with respect to 11, U2, U3, We 


can write 
A= Ao + A uy + Aolts + A 33 
+ (Asis + 2A yottitt2 + Agottg +°+- + A 33443), (3) 
and for B, C, D, correspondingly. Inserting (3) in (A) and comparing coeffi- 
cients leads to the first-order equation: 
AgpDo — BoCo = 1; (4) 
the third-order equations 
AD, + A Do = BC + BiCo, 
ApDe of ADo = BC - BL, (Sa) 
AogD3 fe AzsDpo = BoC3 + B3Co; 


and 
Ao(C2 — Di) — Co(A2 — Bi) = 0, 
Ao(Cs — Ds) + BelCe — Ds) — ColAa — Bs) — DolAy— Br) = 0, (5b) 
Bo(Cs — Ds) — Do(As — Bz) = 0; 


and finally the fifth-order equations: 





* If we differentiate the finite equation above with respect to 1, uz, “3 and subtract from 
each of the equations (7) of the former paper,' equation (A) above results. 
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AD + DoAn — BCu — CoBu = 2(BiCi — AD), 
ApDiz + DoAie — BoCye2 — CoBiz = BiC2 + BL1 — AW2 — AD, 
ApDiz + DoAis — BoCis — CoBis = BiC3 + BsC;) — ADs — AsDi, 
AoDoo + DoA22 — BoCo2 — CoBoe = 2(BxXC2 — AD2), 
ApDoz + DoAo3 — BoCo3 — CoBos = BCs + BsxC2 — A2D3 — AsD2, 
AoD33 + DoAs3 — BoCs3 — CoBss = 2(BsCs — AsDs), 


(6a) 


and 
Ao(Co1—D11) —Co(A 21 — Bus) =Ci(A 2— Bi) —A 1(C2— Ds) — 2(A1C2- AL), 
A o(Co2—D12) —Co(A 22— Bis) =C2(A2— Bi) —A2(C2—D1) — (A D2—A aD) 
—(B,C.—BX)), 
Ao(Co3—D13) —Co(A 23— Bis) =C3(A2— Bi) —As(C2—D1) — 2(BiD2— BaD), 
Ao(C31— Doi) + Bo(C21— Dir) —Co(A 31 — Bas) —Do(Au— By) (6b) 
= — By(C.—D;) — A (C3 — D2) + Di(A2— Bi) +C1(A 3— Bo) — 2(A1C3—A C1), 
A o(C32—De2) + Bo(C22— D2) —Co(A 22— Box) — Do(A22— Biz) 
= — Bo(C2—D;) —Ao(C3—D2) + D2(A2— By) +C2(As— Bs) 
—(A\D3—A3D;) —(BiC3— BxC)), 
Ao(Ca3— Dos) + Bo(C23—D13) —Co(A 23— Bos) —Do(A 23— Bis) 
= — Bs(Co—D;) —A3(C3—D2) +D3(A2— Bi) +C3(A 3 — Bz) — 2(BD3— BsD), 
Bo(Cs,—D21) —Do(Au1— Bu) = D1(As— Bz) — B1(Cs—D2) —2(AL3— AL), 
Bo(C32—Do2) —Do(A32— Box) = D2(A3— Ba) — Bo(Cs— D2) 
—(A2D3—A3D2) —(BX3— BxC2), 
Bo(C33—De3) — Do(As3— Bos) = D3(A3— Be) — Bs(Cs— D2) — 2(B2D3— B;De). 
Moreover, the (6+9) equations (6) are not independent; they are con- 
nected by the identity 
[4 o(Cs2—De2) + Bo(C22—D12) —Co(A s2— Baz) — Do(A 22— Bi») | 
4+ [A pDoo+DoA 29— BoC22—CoB22] — [A Dis +DoA 13— BoCis—CoBis] (7) 
— [Bo(Cn—Da1) —Do(As1— Ba) ] — [A o(C2s— Dis) —Co(A 23— Bis) ] =0. 
2. Gaussian optics. Let us consider first the rays in the neighborhood of 


the axis. Let us assume that 1, “2, us are so small that we can assume func- 
tions A, B, C, D to be equal to their constant members: 


x’ = Agu + Bot,  & = Cox + Dok, 
y = Aoy+ Bon, =’ = Coy + Don; 


(8) 


where ApD»— BoCo=1. 








72 M. HERZBERGER [Vol. I, No. 1 


The evaluation of these equations and the investigation of the geometrical 
meaning of the coefficients form the content of Gaussian optics. 

We shall not go into great detail here, but refer the reader to the discus- 
sion in the Journal of the Optical Society of America. 

Equations (8) can be inverted, and we obtain then 


x = Dox’ — Bot’, E = — Cox’ + Agt’, 
y = Dov’ — Bar’, n = — Coy’ + Aon’. 


(9) 


Let us investigate what happens if one of the coefficients vanishes. 

Ay=0 means that for £=)=0, x’=y’=0, which means that the bundle 
of rays parallel to the axis converges to the image origin. We say that the 
image origin is at the focal point of the system. 

By=0 means that for x =y=0, x’=y’=0. The rays through the object 
origin meet at the image origin. We say then that object and image origin are 
optically conjugate. 

Cy=0 means that £ =7 =0 implies £’ =n’ =0, or, a bundle of rays entering 
the system parallel to the axis emerges parallel to the axis. The system is a 
telescopic system. 

D,=0 means that x = y =0 implies £’=7’ =0. The rays through the object 
origin emerge parallel to the axis. The object origin is the object-side (front) 
focal point. 

3. Image error functions. Let us for finite rays* project image point and 
direction back into the object space, according to Gaussian optics. That 
means we form equations (9) for our finite rays. The ensuing expressions may 
be called the equivalent object coordinates £, 9, —, 7. We have from (9) and (1) 


= = Dox’ — Byt! = (DoA — BoC)x + (DoB — BoD)E = ax + BE, 


? (10) 
E = - Cox’ + A oé’ =(-— CoA + A (C)x + = CoB + AD)é = €X + dé; 
and analogously, 
y = ay + bn, 
n = cy + dn. 
a, 6, c,d are with A, B, C, D functions of 1, ue, u3, and we have 
ao => do => i. 
(11) 


bo = Co = 0, 


the values do, bo, co, do, being the limits of a, 6, c, d for u;=0. If Gaussian 
optics were correct, we would have equation (11) for all values of u;, that is 
for finite aperture and finite field. The deviation from its constant term as a 


2 M. Herzberger, On the fundamental optical invariant, the optical tetrality principle, and on 
the new development of Gaussian optics based on this law, J. Opt. Soc. Amer. 25, 295-304 (1935). 
* The expression finite is used in distinction from paraxial rays, rays near the axis. 


1943] A DIRECT IMAGE ERROR THEORY 73 


function of aperture and field is therefore a measure of the image errors. We 
call a, b, c, d the error functions. 
In the nomenclature of matrix algebra we can express these equations as 


x’ x xo & 
(MC) Ge) = ee) 
E E fo E 
y’ y yo y 
(’)-W(”) (%)=m(?) 
n n no n 
xé, yd, &, nd would be the coordinates of the image ray if Gaussian optics 
were valid. Equations (12) combine equations (1) and (8), M being the matrix 


A B . . . P ee Ao Bo 
S na and M, being the matrix (4' 4 


We have then 
x ~i x x 
is) = My u( ) m( ), 
E 3 g 


s=eML MN. Mo Men. (13) 


follows: 


Let 


(12) 


II 


where 


From (13) it is obvious that the determinant of m is equal to unity. 
Therefore, 
ad — bc = 1. (14) 
The reader can verify for himself that a, 5, c, d fulfill equations (A) and 
therefore equations (5) and (6), which simplify considerably, owing to the 
fact that (11) is fulfilled. 
Equations (13) can be written explicitly 


a= DA — BC, A = Aya + Boc, 
b = DoB — BD, B = Aob + Bod, (15) 
c= -CA+ACL, C = Coa + Doc, 
d= —C,oB+ AWD, D = Cob + Dod. 


Differentiation and substitution in (A) prove our statement. 
4. Third-order theory. The third-order image errors are usually called 
Seidel errors. 
From our point of view, we obtain the image errors by inserting (11) into 
(5). Abbreviating (0a/0ux)u,<0 by a 
a4, +d,=0 Co = dy 
a2 + d,= 0 = by (16) 
a3 + ds; a3 = bo. 


I 
° 
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Equations (16) lead to the conclusion that only six of these twelve coeff- 
cients are independent. Equations (16) are identically fulfilled by selecting 
six parameters k,, with permutable indices such that 


a; = ko bi = kai a=—ky d, = — kay 
de = Ree bs = kao Co = — kis dy = — koe (17) 
a3 = kos bs = kss C3 = — kis dz = — kos. 


Geometrical investigation (which we omit) would show that (if object and 
image planes are optically conjugated), ks3 may be interpreted as the coefhi- 
cient of spherical aberration for the object origin; k23 as the coma coefficient; 
koe and ki; as coefficients of the field errors; and ki as the coefficient of the 
distortion for an object at the origin and an infinite entrance pupil. 

On the other hand, ky, may be considered as the coefficient of spherical 
aberration for an infinite object; kiz as coma coefficient; kis and kz as field 
errors; and kos as the coefficient of distortion for an infinite object and the 
entrance pupil at the object origin. 

The connections between these errors are well-known laws of the Seidel 
theory. 

The method developed here differs from the usual methods in that, first, 
we do not assume the coordinate origins to be in conjugated planes, and, sec- 
ond, we do not restrict ourselves to the consideration of the deviation of the 
object point, but investigate at the same time the deviation of the direction 
of the ray. Equations (10) give, within the limits of our Seidel region, the 
following equations: 

% — x = (hoytty + kootte + kosus)x + (Raitti + kaos + kasus)é, 


& —E = (hiyui + Rims + Risus)x + (Rois + koous + hogtts)é; 


and y and 7 analogously. 

We recommend a detailed study of these equations and their derivatives 
with respect to x and & for meridian rays (y=n=0), especially in the case 
where our origins are not conjugated. 

5. Fifth-order aberrations. For the fifth-order aberrations we find from 
(6a) and (6b) the following fourteen independent equations between the 
twenty-four coefficients a,,, etc. 

Making use of equations (11) and (17) we find that 
do9+do2=2 bum ki2ko3), 


( 
do3+d23= 2 ko2ko3— kiskos— kiokss, 
(ke 


(18) 


Qy+dy=2(ki2— kirksa), 
432+d12=2hi2ko2— kiokis— Riskos, 





Ona das hilig beg ~ Bastien, 
6 — dy = 2hiet+ kirkis— 3kiskes, 
a a) ar a oy oe Oe 
6sa~dag™ Baa} Bishes— 2B rehee, 
6514-11 ™ 3( Biskis— hurkes). 


a, er ee (19) 
boo— do3 = 2kiskes— kiokss— Reokes, 
bos— a33= 2kost+ kisks3— 3ko2ks:, 


C33 +b13= 3(Riskes— kiokss). 
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Here again, we can express the twenty-four quantities in terms of the 
third-order coefficients and nine parameters k,... 

6. The single sphere and the plane. In our previous paper we were able 
to calculate the functions A, B, C, D for the case of a plane and a sphere. 

In the case of a plane, we put object and image origins at the point where 
the axis intersects the plane and found that 


/ , 

x = 2, y = y; 

. , (20) 
f= §&, 7 =; 


or A=D=1, B=C=0. In this case we have a=d=1, }=c=0, and all the 
image errors vanish. 

In the case of a spherical surface, we put the object and image origins at 
the center, and found that 














zs’ =A x, y’ =A yy 
(21) 
= Cx+Di 1 =Cy+Dn; 
where 
1( / 2n®u,—12 2. aS 
C= {4/1 ie domes 2 sf 
r \ r* a 
ie (22) 


2n°U;— U2 {{ ,, 22*u,—% , 2n?u,—U 2 
n2D = +H/ n’'2?— - a n*— - — }—C Uo, 
9 9 9 
r? r- P= 


A= 





D- 


If we develop A, B, C as functions of 4, ue, v3, we obtain the Seidel and 
fifth-order coefficients. Taking care of (5) and observing that 


n 
Ap =-— ’ B, = 0, 
n! 
(23) 
n'’—n n’ 
Co = ——— Do=—> 
r nN 


we finally find the image-error coefficients: 


n? {1 1 \? 1 1 1 
a= ( - P gah eee Seen e's, a3=0, 
r? \n’ n r \n' 


=( 1 4 4 5) n? { 1 1 
qe f — 4), 
r>\n% nn n'n? n3 P 
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and the fifth-order coefficients: 


n* { 1 L\77 3 2 
| doaiaeee ‘Rear tacteees eae <7 
r4§ \n n n* nN 


n* {1 1 2 3 
oe, Pagh TT cee os —+- 
r nN nN nN n'n 





11= b12= doo = b13= bo3= b33= 0, 


~~ ( 1 1 
ri \n!' nn 


645™ 
n*/ 1 1\2/ 2 
c= se reg on aa 
rr? \n n n'> n'n 
n* {1 1 Pe 
SS coin ee Steere 
r?> \n n n' n'n 
~( 1 1 
ig ——{ ———— ] , 
rié\n'2? xn? 


)( 1 3 
nit '3 
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we 
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a13=0, 
do3 >= Q, 
a33> Q, 
1 
a ’ C13=90, 
n4 
(24b) 
Co3=0, 
C33=0, 
d)3=0, 
do3= 0, 
d33=0, 


equations which fulfill all the conditions of equations (6). 


The nonvanishing seventh-order coefficients for one surface 


—( 1 1 )( 5 6 5 
gee ee = ~— fi, 
76 \n!’ nn n'4* n'y mn’ n4 


n* 


(- 1 
"gi Ny! on 


Te ] 1 
ma —T 
i n n nn 
am®/1 1\/5 17 206 33 
—_— ae / : ; 16 = 15, ‘ 14 —_ 13, 
r n nN n nn mn" nn 


10 6 


1'2n? = 'n? 


19 


3 
? +=), 
7 n° 


would be: 
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2 17 5 
5 ant wad tt)” 
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n® {1 1\?/ 8 19 25 19 8 
ee & be =] Cora a3) 
3 7 11 7 3 
(= 7 nn ne sata) 
3n* {1 Ly I 
_— (>-- ) n'n 
3n*® / 1 L477 i 2 2 Z 1 
a G oe, ee ee. 
’ 1 1 1 1 
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oy | 2 1 
~. er ae ’ 
n’? n'n n? 





(25) 
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THE IMPEDANCE OF A TRANSVERSE WIRE IN A 
RECTANGULAR WAVE GUIDE* 


BY 
S. A. SCHELKUNOFF 
Bell Telephone Laboratories 
The purpose of this paper is to derive approximate formulae for the im- 
pedance of a transverse wire carrying uniform current (Fig. 1). 
The total impedance Z to the current through the wire may be defined as 


V 
Z=—) (1) 


where V is the applied voltage and J is the electric current in the wire. The 
total electromotive force V is the sum! of the internal electromotive force V; 
and the external electromotive force V, 
V=V;+ V.. (2) 
Correspondingly, we have an internal impedance Z; of the wire and the ex- 
ternal impedance Z,. By (1) these two impedances are in series with each 
other 
Z=LZi+2Z.. (3) 
If the guide is infinitely long on both sides of the wire, the external im- 
pedance (above the absolute cut-off frequency) is complex 
Ze= Re t+ 1X. (4) 
The resistance term represents radiation of energy into the guide. If the 
frequency is within the principal frequency range and if K is the character- 
istic impedance of the guide to the dominant wave, as seen from the wire,’ 
then evidently 
R, = 4K. (5) 
We shall now calculate the impedance of the wire on the assumption that 
its radius is small. The current in the wire will generate transverse electric 
waves in which the field is independent of the y-coordinate. The general form 
of the field (for z>0) is 





* Received February 8, 1943. 

1 For an explanation, see S. A. Schelkunoff, Electromagnetic Waves, D. Van Nostrand Com- 
pany, Inc., 1943. 

2 And not from a plane current sheet generating a pure dominant wave. 
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~- lrx 
H A2z,%) = be H, sin — eT iz, 
l=1 a 
(6) 
ae _ lex 
— mz K,f14, sin — eT vz, 


l=1 a 


E, (x, 2) 


where I’; and K;, are respectively the propagation constant and the specific 
impedance of a typical TE;,o-wave 
l?x? Ax? : twp - 
dead | arr ale the os hie il (7) 
a“ ? I 1 








The propagation constant of the dominant 7, ,o-wave is 


wr) Ar? Qri mm 
fiw ghee hee (8) 
a? 2 r 4a? 


The dominant wavelength range extends from \,;= 2a to \2=a, Az being the 
cut-off wavelength of TE2,o-wave. If a<\<2a, all the propagation con- 
stants of secondary waves are real 
ee ae (9) 
,=— -—, >t. 
a 17)? 


For the specific impedances we obtain 


h2 \ 1/2 7" 
mae age 
4a’ € 


. iwpa 4a*\—"/? _ 2a 4a*\—1/2 
K, -{1— ) = in —{ 1 — — ‘ 
lx 1?)? Ir 12)? 


The external electromotive force V, necessary to support current J 





(10) 








through a thin wire of radius 7 is 
V.= — dbE,(d,r), (11) 


where d is the distance shown in Fig. 1. This equation presupposes that the 


| | 


= am ee 


‘tt 


Fic. 1 Fic. 2 
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current is distributed uniformly around the wire. As the radius of the wire 
increases, the current distribution gradually begins to depart from uni- 
formity. From (1), (6), and (11), we have 
V ‘3 lxd 
Z.= —=—)> KA; sin — tr, (12) 


I f jul a 





The next step is to calculate the coefficients H;. We shall assume that the 
wire is so thin that the field outside the wire could be regarded as nearly equal 
to that of an infinitely thin electric current filament along the axis of the wire. 
For an infinitely thin filament, we have 


 &£fory . tes 
H, = lim —f — sin —dx as r— 0. (13) 
advq 4 a 


Integrating and passing to the limit, we obtain 





I . Ind 
H, = — sin —- (14) 
a a 
Substituting (14) in (12), we have 
= lad 
oe™ — >> Ki sin? — eT; (15) 
2 ji=1 a 
and, therefore, 
b wd 
R,. = — K,sin?—» 
a a 
2b we 
K = 2R. = — Ki sin? — » (16) 
a a 
's lxrd 
iX, = — )> Ki, sin? — er. 
a l=? a 
Substituting from (10), we obtain 
: 2b . «ad A? \71/2 
K = y—sin?’—{ 1 —- — 
a a 4a? 
(17) 
. 2b > 1 (: i R , ind . 
X,=1— — = sin? —— e-Tir, 
r l=2 l J?) a 


Hence, the ratio of the external reactance of the wire to the characteristic 
impedance of the guide (as seen from the wire) is 


X. a ON a 2, 1 4a\-/2 Id 
—_— = — 1—— csc?— > —( 1 —- — sin? —— eT, (18) 
K r 4a? a ime / P)? a 
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It is evident that the total inductance of the wire is a series combination 
of inductances associated with the individual secondary TE waves, generated 
by the current in the wire. 

For the internal impedance of the wire, we have® 

bIo(o;:7r) 


2rr1,(o;7) 


Sy iwi (gi + iwe;), n= —_—__——— 
gi + Wwe; 


This is a general expression applicable to dielectric wires as well as to metal 
wires. In the case of metal wires, we let e;=0. Usually, the radii of metal 
wires will be sufficiently large to make the modified Bessel functions in (19) 
nearly equal so that approximately 





Zi = Ni 


(19) 





bn; b ip; b has 
Bp oe ti —/— = 4/7 (1+ 3). (20) 
2ar ar 2: 2ar gi 


If r/a is small, the series (18) converges slowly. The difficulty may be 
obviated with the aid of the following device. Let 


_ = p>» uy (21) 


be a slowly converging series; let 


v= dD» (22) 
l 


be a series of terms approximating (21) in such a way that the approximation 
becomes increasingly better as / increases; then 
u= dint Do (um — 0) (23) 
l l 

so that (21) can be regarded as the sum of two series, of which the second is 
more rapidly convergent than the original series. If now the sum of the first 
series in (23) happens to be known, we have succeeded in replacing the 
original slowly convergent by a more rapidly convergent series. 

We shall apply this device to (18). First we rewrite (18) in the following 
form 


lard 
1 — cos—— 





xX. 1 fia rd = a 
—- = — post —1csc?— >> —— or (24) 


K d? a =2 a* 
4 4 l 1 / | 
17)? 


3 See the book mentioned in footnote 1. 
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then we note that as/ tends to infinity, I’; tends to /r/aand V1 —4a?/Pd? tends 
to unity. Hence, a typical term of the v-series will be (1/)[1— 
cos(2/rd/a)|e~'*"/*, and (24) may be expressed as 














2lad 
x ' = q 1 — cos—— 
he /4a? T < a 
— = —4/—-1 csc? | SS rive + r| 
K 4 ? a l=2 t 
(25) 
2lad 
1 — cos 
os a e Tir 
T = ee cir), 
D l ( j 4a? 
ian 
It is known that 
1 1 
} & 3 e—'? cos lq = — a log (1 — 2e-? cos g + e~??); (26) 
i=1 4 
therefore, 
 - 1 /4a? rdf 1 2nd 
“thee ——— 1 cs — E log (1 — 2e-*"!* cos — + e*evl* 
K 4 d? a 2 a 
2rd re 
— log (1 — e-*"/*) — eme(1 — cos ) +T}. (27) 
a 


This can be transformed into 


SS ee: rd[ 1 Sa 2ad 
—$ = — —1 eset =| — log (cosh —_—- cos") 


K iV ? a a a 





1 Tr 2xd m 
— — log 2 — log sinh — — eel —cos—)]+T7]|. (28) 
? da a 


An entirely different expression for Z, can be obtained by the image 
method. Assuming again that 7 is sufficiently small and that the wire is not 
too close to the boundaries of the guide and that consequently there is no 
“proximity effect,” we can immediately obtain 


1 ; = 2 _ 
Ze= = nBb | wetan + 2>> H,(2nga) —  & H,(2nBa + 28d) 


n=1 n=0 
— > Ai(2nba + 28a — 28d) |. (29) 
n=0 


This is a slowly converging series and is useless for direct numerical com- 
putations; on the other hand, it may be useful for other purposes. Thus the 
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difference between the external impedances of two wires of different radii is 
obtained in the following simple form 


Z Are) —Z.(r1) = 1n8b[Ho(Br2) a H4(6r1) | 
” inBb [J o(Bre) - J (Bri) | + 4n8bi[No(Br1) ses No(6r2) J. 


The first term represents the effect‘ of the radius of the wire on the impedance 
of the guide as seen by the wire. The second term represents the difference 
between the external reactances of two wires 


X.A(r2) — X.(r1) = 49Bb[No(Bri1) — No(Bre) |. (31) 


(30) 


This equation can be used for numerical calculations in conjunction with (28). 
The slowly converging part of (29) is the mutual impedance between the wire 
and the wave guide. 

An expression for the mutual impedance between two parallel wires in 
the wave guide can also be obtained by the image method. Thus we have 
(de >d;) 


1 : 
Zin = = 8b | wsiads — Bd:) + D> Ho(2nba + Bd — Bay) 


n=1 


+ Y>Hi(2nBa + Bd; — Bd2) — D> Ho(2nBa + Bas + Bd) 


n n=0 


= 2 
— >> Ho(2nBa + 28a — Bd; — ads). (32) 
n=( 

Next we shall deal briefly with the external impedance of a “split” wire 
(Fig. 2). Let the electromotive intensity at the surface of the wire and the 
current in the wire be 








mn y 
E, = — )> Em cos ——: 
= (33) 
“= mnry 
I= > I, cos : 
m=0 
The complex flow of power is 
». Beene wt 3 + 1s 
Vv al Saclad idame” Eolo + > 2. Emlm |, (34) 
m=1 


where Y, is the external admittance of the split wire and V, the voltage across 
this admittance. Since 


V. = bE, (35) 


4 Which was entirely ignored in the derivation of (28). 
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we obtain from (34) 


= 0 1 = > yee 
oy (36) 
- Eb deat EoE& 
where 
a Pe ate 
Y.. (3/) 


bEm 

The first term in (36) is the external admittance corresponding to a uni- 
form current filament (Fig. 1) and, hence, is equal to the reciprocal of either 
(15) or (29). The second term is the capacitative admittance (assuming that 
\>2b) between the external surfaces of the two portions of the transverse 
wire. The impedance diagram is shown in Fig. 3 where the parallel lines 


LJ 




















Fic. 3 Fic. 4 
represent the wave guide, the inductive reactance X; is the reactive part of 
(29), the capacitative reactance X{’ is the reciprocal of the second term in 
(36) and X; is the reactance looking inward from the surface of the gap in the 
wire. In the case illustrated by Fig. 4 the internal reactance is approximately 
a s 
Ai= ge (38) 
iwemr? 


where r is the radius of the wire and s is the length of the gap. 

The internal reactance of two hollow cylinders as well as the quantity 
X 2’ will be discussed in a separate paper. Here we shall merely derive general 
formulae and show that roughly X2’ is equal to the external capacitative 
reactance between two sections of a split transverse wire placed across infinite 
parallel planes. 

Each component J, cos (mmy/b) of the total current J in (33) originates a 
radial wave of the following type 


I mK 1(¥ mp) mn y 
H 4(p) aes QarK C. —— COS > 
£mTr mY 
seals (39) 
NY ml mK o(¥ mp) mny 
[oe oo 





Bd an tus ee 
v(e) 287K 1(ymr) b 
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where p is the distance from the axis of the current and y, is the radial 
propagation constant of the mth cylindrical wave 


mx? = Arr? mr a 
Ym = —— ea ee SS 1 -— ° (40) 
2 2 b m?)? 





If \> 2b, all the radial propagation constants of order m higher than zero are 
real. This explains why even the nearest image will have but little effect on 
the admittance Y, except when the wire is quite close to the walls of the 
wave guide, or when J is nearly equal to 2b. Even when the wire is close to 
the walls of the guide only the nearest image will have an appreciable effect 
on Y,, unless \ is nearly equal to 20. 

The complete expression for the impedance Z,,=1/ Y» is 


n=1 


NY md 
lm — _ Sie yas | Keyan + 2 +> Ko(2nyma) 
2riBrK 1(y mr) 


as Ko(2nyma + 2ymd) — bbs Ko(2nyma + 2y¥ma + 2ad) |. (41) 
r=0 


n=0 


Equation (29) corresponding to the principal cylindrical wave (m=0), is 
of course, a special case of (41). The propagation constant I» of the principal 
wave, however, is pure imaginary and, hence, distant images have a pro- 


nounced effect on Z,. 
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NOTES 
ON A. C. AITKEN’S METHOD OF INTERPOLATION* 
By WILLY FELLER (Brown University) 


1. A. C. Aitken! has recently devised a method of practical interpolation 
which is particularly well adapted for computing machines; neither differ- 
ences nor tables of interpolation coefficients are used, and the necessary opera- 
tions are most easily performed on modern computing machines. Moreover, 
the degree of the interpolating polynomial decreases by two at each stage, 
which minimizes the amount of necessary work. Recent experience has again 
confirmed that the method is extremely convenient and timesaving. It would 
nevertheless seem that the method is not sufficiently known, and we propose 
therefore to give a brief outline. Our proofs seem simpler than the two proofs 
given by Aitken,! or that given by Lidstone.* At the same time we shall be 
led to a procedure which works for an odd number of data as well as for an 
even number (originally the method appeared to work for an even number 
only and special computational devices were used to reduce an odd number 
of data). For most practical arrangements of computations we have to refer 
to Aitken':? and Lidstone.® 

2. Linear cross-means. The full power of the method appears only with 
the use of quadratic cross-means, but these are in turn based on linear cross- 
means. Moreover, with completely unsymmetrical data only linear cross- 
means can be used. 

Let it be required to compute the value f(£) of a polynomial of mth degree, 


f(x), given f, =f(x.) for R=0, - - - , nm. We note that 
fo = | . 
f(x) =| ~ “s +(e ~ a) (1) 
f(z) z-€ 


is a polynomial of degree n—1, and that f‘(£) =f(£). Hence we are required 
to find f((&) knowing 


( | f x0 = é | ; 
rg = F |= (x — Xo) (2) 
tk Fo E 
for k=1,---, . Thus the problem has been reduced from m to n—1. In 


like manner the problem is further reduced to n—2 and so on. 
All the computer has to do is to write in a column the “parts” x,—£&, and 


* Received Dec. 18, 1942. 

1A. C. Aitken: On interpolation by iteration of proportional parts, without the use of differ- 
ences. Proc. Edinburgh Math. Soc. Ser. 2, 3, 56-76 (1932). 

2 A.C. Aitken: Studies in practical mathematics III: The application of quadratic extrapo- 
lation to the evaluation of derivatives, and to inverse interpolation Proc Roy Soc. Edinburgh, 
Vol. 58, pp. 161-175, 1938. 

*G. J. Lidstone: Notes on interpolation. J. Inst. Actuar., Vol. 68, pp. 267-296, 1938. 
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in an adjacent column the corresponding values f;. Using (2), new columns 
are successively added to the right, the number of rows decreasing by one 
each time. The “parts” remain obviously the same throughout the computa- 
tion. The determinant in (2) is easily computed, and the result appears in 
the main dials ready for division without clearing. Moreover, on most ma- 
chines, the divisor x, —xo = (x, — £) — (xo — £) will automatically appear on the 
secondary dials (provided the main keyboard has been used for the factors f;). 
Actually in most cases the divisor x,—xo will be a small integer. It should 
also be noted that as the computation proceeds the entries will tend to agree 
in an ever increasing number of their more important digits. These, of course, 
will not be copied down; this reduction of digits of f, makes it in turn possible 
to drop some last digits of the “parts.” 

3. Quadratic cross-means. For these it is necessary that the given data 
be placed symmetrically with respect to some point x =a. Denote, then, two 
symmetrically placed points by x, and x, (R=1,---, m; x,.-a@=a—x_,). 
The point x» =a is included among the data only if m= 2m. Consider 


| f(2a — x) 2a—x—€| 


» = +is~— 3 
(x) | f(x) o<% | (x a) (3) 
and 
| — f(2a — aa~2a-—§ 1 
y(x) = | HK ' *) : ‘ | 4- 2(é —- a). (4) 
| f(x) on 


Obviously $(x) and W(x) are even functions of x—a, and hence polynomials 
in t=(x—a)*. Moreover, $(&) =W(&) =f(&). 

(a) If m=2m-—1, the problem is reduced to finding the value for ¢ = ( —a)? 
of the polynomial of (m—1)th degree ¢(a++V/2) given its values 


| fe *-e — € | : 
Pan : | + (x, — x2) (5) 
Si te — & | 
for t=(x,—a)*, R=1, ---, m. Thus a simple application of (5) will reduce 


the number of data from 2m to m. From here we proceed as before using 
linear cross-means. It should be noticed that the “parts” now to be used are 
(x, —a)*—(£—a)? = — (x,—£)(x_.—£), that is to say the product of the parts 
already used for (5). This invariant property dispenses of the necessity to 
label the panels. In most practical cases the new “parts” will differ only by 
integers or multiples of 1/2. 
(b) If »=2m, we compute the values 
Po ewe Se (6) 
Se te—& 

for t=(x.—a)?, k=0, - - - , m, and proceed as before. Since here the denomi- 
nator is the same for all & the division may be deferred to the final result. 
This is a slight simplication. 
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A NEW DERIVATION OF MUNK’S FORMULAE* 
By W. C. RANDELS! (University of Oklahoma) 


Recently M. A. Biot? has applied the method of the acceleration po- 
tential to some problems of two-dimensional airfoil theory. In this paper this 
method will be used in order to obtain a short proof of Munk’s formulae?’ 
for the lift and moment of a thin airfoil. 

As usual in the theory of thin wing sections we replace the airfoil by its 
mean camber line supposed to deviate but little from the chord. Studying the 
plane irrotational flow of an incompressible fluid around this indefinitely thin 
airfoil, we take its chord as the x-axis of a system of rectangular coordinates 
x, y, ascribing to the leading and trailing edge the abscissae —1 and +1 re- 
spectively. Denoting by V the velocity at infinity and by a@ the angle of at- 
tack, supposed to be small, we write the x- and y-components of the velocity 
vector was V+u and aV-+v respectively, where u, v and aV will be small 
as compared with V. We denote the pressure by p and the density by p. 


Then, by Bernoulli’s equation 


.~) 





w? + p = po = const. 


bo 


Neglecting quantities of the second order, we have 

p(V? + 2uV) = — 2(p — fo). (1) 
The quantity ®= —1/p(p— po) is called the acceleration potential, since the 
acceleration equals grad ®. 

It is known that V-+u—1i(aV+v) is an analytic function of z=x+1y. 
Since V is a constant (}V?+uV)—i(aV*+vV) is also an analytic function 
of z. The functions ®=(3V?+uV), V=—(aV?+vV) are thus seen to be con- 
jugate harmonic functions. 

Let the mean camber line of the airfoil be given by the equation y =-c(x), 
(—1<x<1; c(1)=c(—1)=0). The condition that this be part of a stream 
line furnishes the condition 

aV+u : : 
— = '(x) along y = c(x), (-1S x8 1). 


V+u4 





* Received Oct. 28, 1942. 

1 This note has been prepared at the suggestion of Professor W. Prager while the author 
was a fellow under the Program of Advanced Instruction and Research in Mechanics at Brown 
University. The author is indebted to Dr. L. Bers for valuable advice. 

2M. A. Biot, Some simplified methods in airfoil theory, Journal of the Aeronautical Sci- 


ences 9, 185-190, (1942). 
3M. M. Munk, General theory of wing sections, N.A.C.A. Techn. Rep. No. 142 (1922). 
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Neglecting quantities which are small of a higher order than the first, we 


obtain 
VW = — V*c'(x) along y = c(x), (-1S5 x81). (2) 
Since v vanishes at infinity we have 
W(o2) = — aV?. 


As the mean camber line deviates but little from the segment —1<x<1 
of the x-axis, we will not commit an appreciable error by fulfilling the condi- 
tion (2) along this segment rather than along the mean camber line. We set 


v= -—aV?+ U4 ¥2 
where 


VY, = aV? and ¥, = — V*c'(x) along — 15 x $1, y=0 


and 
V (2) = W2( 0) as (, 


VW, and the conjugate harmonic function ®, have been determined by Biot. 
From ®, the lift distribution due to the angle of attack can be obtained. In 
the following we shall set a=0 and thus obtain the lift distribution due to 
the curvature of the mean camber line. Within the framework of our linear 
theory these two influences are additive. 

In order to solve the boundary value problem for Wz we map the exterior 
of the segment of the real axis between z= —1 and z=-+1 onto the exterior 
of the unit circle in the ¢ plane by the conformal transformation 


z= 1/2(¢ + 1/f). 


The line segment (—1, 1) then is transformed into the circumference of the 
unit circle and we have x =cos @ (Fig. 1). Since a conformal transformation 
takes a harmonic function into a harmonic function, our problem becomes 
that of finding a harmonic function having the values — V%c’(cos @) on the 
unit circle. If we assume V to be regular on the boundary, the solution is 
given by the Poisson integral but the resulting function will not vanish at 
infinity unless 


2 


1 
ayo = — c’(cos 0)dé = 0. 
2r 0 
Therefore, to satisfy the condition Y( «) =0 we introduce a singularity cor- 
responding to a source-sink doublet at the leading edge.‘ With the notations 
of Fig. 1, we obtain 
‘It is natural to assume a singularity at the leading edge since our assumption about u 


and v being small does not hold here. 








90 NOTES [Vol. I, No. 1 
cos 6; 1 r? — 1 


V(r, 0) = — 2a.V?- -—vif c'(cosr) — a dr. 
( . ry 2r " 0 le" ee Per gee 


This function clearly vanishes at infinity and will satisfy the other boundary 
condition because 








cos 4; 1 —— 
= = on the unit circle. 





rT) 








3 
G 
! 
! 








Sex) 


x 








Fic. 1. 


The Poisson integral used above is only legitimate if 


Qe : 1 | c’(x) | 
| c’(cos r)| dr = 2 penn EE 
0 1 [1 — x?}1? 


This implies a condition on the rapidity with which c(x) tends to zero as x 
tends to +1. 

The values of the conjugate function ®(r, 6) on the boundary of the unit 
circle are given by the formula® 


sin 0; 1 ail 7— 96 
/ " 
— c’(cos r) cot ——— dr, 
T1 2r 0 2 


$(1, 0) = — 2aV? 





where /” denotes the Cauchy principal value. The total lift L will be given by 





5 J. D. Tamarkin, “Theory of Fourier series,” Brown University, 1933, p. 110. 
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= 


2r 
_= rf ®(1, 0) sin 6d0 
0 


2? sin 6, . 
=— 2aypv? f sin 6d0 
0 


r 


V2 Qn : * 2n a r) 
> — f sin oan f c’(cos r) cot ——— dr. 
2 0 0 2 


T 





It is easy to calculate: 


ad sin 6; . 
f ——— sin 6d0 = r. 
0 


rT 


The second integral is evaluated by making a formal interchange of the order 
of integration. This interchange can be easily justified. We then have: 


1 2r *2r T= 6 
f sin oa [ c’(cos r) cot -dr 
2r 0 0 2 


2e 1 * 297 ? T —6 
= f c’(cos r)dr —f sin 6 cot dé 
0 2a 0 2 


and since the function conjugate to sin 6 is —cos @ 


1 + 3e (r — 6) 
-f sin 6 cot ————— d@ = — cos rT. 
0 2 


2a 








Using this together with the definition of a>) we obtain the lift 


” 


owl 


2 vf *(x) Ad. 
= =— Cc s —_—_—_—_— 
p - . 1 ae x? ]2 


2ev*4 (2) ee: fw ae d 
—2 C(x ——————— _ c( x) ————-— ax 
p [1 re x? ]i/2 - - [1 a ac? ]3/2 


- 1 dx 
= 2h fi a-ni-e)" 


We have assumed that c(x) is such that limz.: (c(x)(1—x)-"? =0. The under- 
lined expression is Munk’s formula for the total lift, due to the curvature of 


the wing. 


® 2r 
c'(cos r)dr +f c'(cos r) cos rar] 
0 
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A similar procedure furnishes the moment M of the lift. It is given by 
M = rf (1, 6) cos @ sin 6d0 
0 
: 2" sin 0; i 
= — 2pV 7a —— cos @ sin 6d0 
0 r} 
pV? Qn * 29 ¢-£ 
+ — cos 6 sin 6d0 c’(cos r) cot dr. 
2a 0 0 2 
Then 
2" sin 0, - 7 
cos # sin @d@ = — — 
0 ri Z 
and 
1 2¢ ’ * 2r ‘T= i] 
— cos @ sin 6d0 c’(cos r) cot ——— dr 
2r 0 0 Zz 
1 2r 
=—-— c’(cos r) cos 2rdr 
? 
= 0 
so that 


1 Qe 1 2r 
M=-— ov] - —f c'(cos r)dr + —{ c’(cos r) cos 2rdr 
y 0 2 0 


) vf " ) a. 1 
= — 2pV? c' (x) ————_—_ dx 
OTe 


, vf (x) xdx 
= Zz . C(x ———_—_—_—_——_——— — 
= | {1 = we? }t/2 


which is Munk’s formula for the moment, due to the curvature of the wing.® 


6 After the manuscript of this paper had been completed (August 1942), H. J. Stewart 
has published an analysis proceeding along similar lines: A simplified two-dimensional theory 
of thin airfoils, Journal of the Aeronautical Sciences 9, 452-456 (Oct. 1942). 
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CONCERNING THE ACCELERATION POTENTIAL* 
By LIPMAN BERSt (Brown University) 


The following lines aim at indicating the possibility of a more rigorous 
approach to Prandtl’s method of the acceleration potential for two dimen- 
sional flow.! 

We consider a steady incompressible potential flow past an airfoil of in- 
finite span. We assume the profile, P, to be given by 


z=x+iy=Z(s), O<s&K sz, (P) 


where the sense of increase of the arc length s corresponds to the counter- 
clockwise direction and the sharp trailing edge, 7, isgiven by T= Z(0) =Z(sr). 
The position of the stagnation point, S, near the nose of the airfoil shall be 
given by S=Z(ss). We also set 


dZ F 
— mn 2 els) | 


8(s) being a continuous function of s and such that on the upper bank of the 
wing near 7, —7/2S8 87/2. 
We denote by u and v the velocity components in the x and y directions 
respectively and assume that 
a= U>0O, v = 0, at infinity. 
Then u—iv is an analytic function of z=x++1y and so is 
© + iv = log (u — iv). 


At S the function ®+7¥ possesses a singularity. (There also is a singularity 
at 7, unless the angle there is 0.) ®+2W may be determined as a solution of 
the following boundary value problem: 


A. To determine a one-valued analytic function ®+2¥ defined on the region 
exterior to P and satisfying on P the boundary condition 


* Received Jan. 22, 1943. 

| This note has been written at the suggestion of Professor W. Prager. The author is 
indebted to Professor K. O. Friedrichs for criticism. 

1 Cf, L. Prandtl, Beitrag zur Theorie der tragenden Flache, Zeitschrift f. ang. Math. u. Mech. 
16, 360-361 (1936); Uber neuere Arbeiten zur Theorie der tragenden Fliche, Proc. Intern. Congress 
for Appl. Mech., Cambridge, U.S.A., (1938), pp. 478-482, (See also: N.A.C.A. Technical 
Memorandum No. 962); M. A. Biot, Some Simplified Methods in Airfoil Theory, Journ. of the 
Aeronaut. Sci., 9, 185-190 (1942); W. C. Randels, A New Derivation of Munk's Formulae, 
this Quarterly 1, 88 (1943). 
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_ a B(s) for OSs Ss (1) 
x—B(s) for sss Sr 
as well as the condition 
& = log U, Y=0, at z= o, (2) 


Equation (1) expresses that P is a streamline of the flow and takes care of 
the Kutta-Joukowsky condition (no flow around 7). The unknown position of 
S is uniquely determined by (2). For instance, let z=f(¢) map [¢| >1 con- 
formally into the exterior of P, taking {= © into z= © and {=1 into z=T. 
Set f(e#®) =Z[o(6)], S=f(e'"). Then the condition (2) may be written in the 


form 


f ¥izlo@ las = 0), 


0 


In view of (1) we obtain 
1 2 
r=2r- —f 8 [o(8) |dé. (3) 
us 0 


From ® we can calculate the pressure p. In fact, we have by Bernoulli's 
equation 
p + 3p(u? + v?) = pu. + 2pU* = po, (4) 
where p is the density, p, the pressure at infinity and pf» the stagnation 
pressure. Since ®=log |u—iv|, 
Pp, " 
p = Po — — e, (5) 
If the wing is infinitely thin, say given by 


S2xs1, (P;) 


IIA 


2 = &, y = Y(z), —1 
the boundary value problem A takes the form: 


B. To determine a one-valued analytic function ®+1V defined on the region 
exterior to P, and satisfying on P, the boundary condition 


‘— arc tan Y’(x) on the upper bank of P,, 


lA 
8 
IA 


v= P.. arc tan Y’(x) — w on the lower bank of P,, Xs 


—1 
— arc tan Y’(x) on the lower bank of P;, xs S$ x 21 
(x3+2Y(xs) being the stagnation point) as well as the condition (2). 


If the wing is very slightly curved and very slightly inclined, the above 
rigorous but inconvenient treatment can be simplified as follows. The dis- 
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tance between the leading edge LD = —1+21Y(—1) and the stagnation point S 
is small of second order as compared with the angle of attack. In fact, the 
general character of the flow around P;, will be similar to that of a flow 


around a straight line, say Pe, 
y= — xtana, — cosa SxS cosa. (P2) 
By 
sep t—ote — (6) 
4 f 


the exterior of P2 is mapped into |¢| >} and T is taken into }e~‘*. S is taken 
into —4e-** (this follows for instance from (3)). Therefore 


S = — Leite + e-tie) 
and, since in this case L = —e~‘*, we have for small values of a 
|L—S| ~ 2a’. 


Now, ®+7Y possesses singularities at Z and at S. For small angles of 
attack we may assume that we will make a very slight error if we replace 
these two singularities by a single singularity situated at LZ. In order to 
determine the character of this singularity, we again consider the flow around 
P,. The complex potential, say for U =1, is given by 


1 
errt.* (i sin a) log ¢ (7) 


so that, by (6) and (7), 
dw dw /dz o + je** 
ds dtl dt fo +4e* 





and 

© + iv = log (¢ + 3e**) — log (¢ + 3e7**). 
This is (in the {-plane) the complex potential of a source at —}e~** and a 
sink at —4e**. For small values of a we may approximate this source-sink 


system by a doublet with a vertical axis. 
Thus we replace problem B by 


C. To determine a one-valued analytic function ®+iV defined on the region 
exterior to P;, satisfying on P, the boundary condition 


W = — arc tan Y’(x) (8) 


as well as the condition (2) and possessing at L a singularity which assumes the 
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form of a potential of a doublet with a vertical axis when the exterior of P, is 
mapped conformally into that of a circle (by a transformation f(z) with f’() >0). 


The actual solution of this problem is still difficult. Therefore we make use 
of the fact that P; is closely approximated by the slit 


y=0, —-1<%*31, (Ps) 


and replace the domain of definition of ®+7¥ by the exterior of P3. Then we 
obtain the following boundary value problem: 


D. To determine a one-valued analytic function ®+1V defined on the region 
exterior to P3, satisfying on P3 condition (8) as well as condition (2) and possess- 
ing at —1 a singularity which, in the [-plane determined by 


(t+) 
z=—(¢+—), 
2 ¢ 


assumes the form of a potential of a doublet with a vertical axis. 


This problem can be easily solved. The presence of the singularity en- 
ables us to satisfy both conditions, (2) and (8). 

It remains to show that the method described above is identical with the 
method of the acceleration potential, the latter usually being presented as 
based upon the assumption 


Ap = 0. 


By virtue of our hypotheses p—p., will be very small as compared to pz — po 
(except at the neighborhood of the leading edge), so that disregarding terms 
of higher than first order in (b—p.)/(P2—po) we have 
; po — Pp 
log (po — p) = log (po — pu) + ——— 


0 «2 


and, by (4) and (5), 
p = — pU*@ + const. 
On the basis of the above considerations an estimation of the error (due 


to replacing the actual problem B first by C and then D) seems to be both 
desirable and possible. 
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